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In the framework of the isomonodromy deformation method, we present a constructive pro- 
cedure to obtain the critical behavior of Painleve VI transcendents and solve the connection 
problem. This procedure yields two and one parameter families of solutions, including trigono- 
metric and logarithmic behaviors, and three classes of solutions with Taylor expansion at a 
critical point. 



1 Introduction 

This paper appeared in May 2006. I put it on the archive now, with more that four years of delay, 
for completeness sake. The paper is published in J.Phys.A: Math. Gen. 39 (2006), 11973-12031, 
with some modifications. The sixth Painleve equation is: 
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(PVI). 



(y- 1 ) 2 (y-x) 2 

The generic solution has essential singularities and/or branch points in 0,1, oo. It's behavior at these 
points will be called critical. The other singularities, which depend on the initial conditions, are 
poles. A solution of PVI can be analytically continued to a meromorphic function on the universal 
covering of P 1 \{0,l,oo}. For generic values of the integration constants and of the parameters 
a,/3,7,(5, it cannot be expressed via elementary or classical transcendental functions. For this reason, 
it is called a Painleve transcendent. Solving (PVI) means: i) Determine the critical behavior of the 
transcendents at the critical points x = 0, 1, oo. Such a behavior must depend on two integration 
constants, ii) Solve the connection problem, namely: find the relation between couples of integration 
constants at x = 0, 1, oo. 

We use a matching procedure to study the above two problems. The procedure allows us to 
compute the first leading terms of the critical behavior at a critical point and the associated mon- 
odromy data. This procedure is essentially the isomonodromy deformation method. The reason 
for our terminology is that we make particular use of the matching between local solutions of two 
different reductions of the linear system of ODE, associated to (PVI) by the isomonodromy defor- 
mation theory. This matching allows us to obtain the leading term(s) of the asymptotic behavior of 
a corresponding Painleve transcendent y(x). In this sense, we say that our approach is constructive. 
Namely, we don't assume any behavior of y(x); rather, we obtain it from the matching condition. 
This differs from other authors' approach, who start by assuming a given asymptotics for y[x) and 
then compute the corresponding monodromy data (and so they solve the connection problem). This 
kind of approach was successfully used for some of the Painleve equations and allowed many pro- 
gresses. Our approach is developed to tackle with the cases when we don't know - or we are not able 
to guess - the asymptotic behavior. In the case of (PVI), we may say that most of the solutions are 
known. But for some points in the space of monodromy data, we still don't know the corresponding 
critical behaviors. Our work is motivated by the need to explore these remaining cases. 

Once the local matching is done, we proceed with a global description of the solutions of the 
associate linear system of ODE, in order to compute its monodromy data. These are the monodromy 
data associated to the solution y(x), of which the asymptotic behavior has been obtained by the 
precedent step. Again, this computation is done by a (global) matching, among solutions of the two 
reduced systems and that of the original one. This is the main powerful point of the isomonodromy 
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deformation method. The monodromy data are computed in terms of the coefficients of the linear 
system of ODE, which are elementary functions of the parameters (namely, the integration constants) 
appearing in the leading term of the asymptotic behavior of y(x). The inversion of the formulae 
expressing the monodromy data, gives the leading term of y(x) in term of the monodromy data. 

The procedure can be repeated at the other singularities x = 1, oo. In case of (PVI), x = 0, 1, oo 
are equivalent by symmetry transformations. These facts allow to solve the connection problem 
([17], [8], [9], [11], [5]). 1 

The work of Jimbo [17] is the first on the subject. For generic values of a, (3, 7 S, PVI admits a 
2-parameter class of solutions, with the following critical behaviors:. 

y{x) = ax 1 - (T (l + 0(\x\ e )), x -> 0, (1) 

2/ (x) = l-a( 1 )(l-x) 1 - ff(1) (l + 0(|l-x| e )), x^l, (2) 

l,(a;) = a(°°^' T<oo) (l + 0(|a;|- e )), x 00, (3) 

where e is a small positive number, and crW are complex numbers such that ^ and 
< Iftcr < 1, < 5Rcr (1) < 1, < 5fter (oc) < 1. We remark that x converges to the critical points 
inside a sector with vertex on the corresponding critical point. The connection problem is to finding 
the relation among the three pairs (a, a), (a^\a ). In [17] the problem is solved 

by the isomonodromy deformation method. In particular, the exponents are determined by the 
relations: 

2cos(tto-) = tr(M M x ), 2 cos(tto- (1) ) = tr(MiM x ), 2 cos(7rcr (oc) ) = tx(M M{). 

Here Mq, M x , M\ are monodromy matrices to be introduced below. 

The above class of solutions was enlarged in [26] and [11], to the values a G C, a £ (— 00, 0] U 
[l,+oo) (here we consider x — > 0). When Kcr > 1 or Ker < 0, the critical behavior is like the 
above, but it holds for x — >• in a spiral-shaped domain in the universal covering of a punctured 
neighborhood of x = 0, along a paths joining a point x to x = 0. Along special paths which approach 
the movable poles, these solution may have behavior y(x) ~ sin - ^ 2 ^ lnx + <fi(x,afj, where <fi(x, a) 
is a phase depending on the parameter a. The transformation a M> ±<r + 2N, ]VeZ, leaves the 
identity tr(M M x ) = 2cos(ira) invariant. Its effect on the solutions is studied in [11]. As a result, 
one can reduce to the values < 5Rcr < 1, a ^ 0,1. The reader may find a synthetic description of 
these results in the review paper [12]. 

It is an open problem to determine the critical behavior, say at x = 0, for a = 0, 1. To be more 
precise, the problem is encountered when tr (MiMj) — ±2. These are precisely the points of the 
space of monodromy data mentioned above, in correspondence of which we do not know the critical 
behavior. In addition, certain non-generic values of a, /?, 7, 5 are not yet studied. The matching 
procedure is motivated by the need to explore these unknown cases. 

As a result of the matching procedure, we obtain: 

Rl) A two-parameter family of solutions, of the type found by Jimbo [17]. Besides, we show 
that there are solutions with trigonometric behavior. 

R2) One-parameter families of solutions, including a class of logarithmic solutions. 
Together with the results of [26] and [11], Rl) and R2) will cover all cases tr(MjMj) ^ —2, namely 
it / 1. By symmetry transformations, some of the cases tr(MjMj) = —2 can be obtained from the 
above results (for example, the Chazy solutions [21]). 

R3) The solutions which admit a Taylor expansion at x = 0. 

R4) We compute the corresponding monodromy data. 

1 For reasons of space, we limit ourselves to the computation of monodromy data, without explaining how the 
connection problem is practically solved once the monodromy data are computed, and how the analytic continuation 
is done. We refer the reader to [17], [8], [9], [11], [5]. The behaviors at x = 1 and x = 00, and the dependence of them 
on the monodromy data are deduced from the behavior at x = by symmetry transformations. PVI is invariant for 
the change of variables y(x) = 1 — y(t), x = 1 — t and simultaneous permutation of $o, 9i ■ This means that y(x) solves 
PVI if and only if y(t) solves PVI with permuted parameters and independent variable t. Similarly, PVI is invariant 
for y(x) = l/y(t), x = l/t and simultaneous permutation of 9^, 9(). It is invariant for y(x) = (y(t) — t)/(l — t), 
x = t/(t — 1) and simultaneous permutation of 80, 6 X . By composing the third, first and again third symmetries, we 
get y(x) = y(t)/t, t = l/x with the permutation of d\, 8 X - Therefore, the critical points 0, 1, 00 are equivalent. 
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In virtue of the symmetries of (PVI) (birational transformations of (x, y(x))), it can be shown that 
the solutions with Taylor expansion at x = 0, obtained by the matching procedure, are the repre- 
sentatives of three equivalent classes, which include all the solutions admitting a Taylor expansion 
at a critical point. If we define a through the relation ti(M M x ) = 2cos(7rcr), the representatives of 
three equivalent classes correspond to values a — 0, a — ± / =F 9 oo) and a = 1. 

A further step in the study of PVI, is the problem of the systematic classification of all the 
solutions of (PVI) in terms of the monodromy data of the associated linear system. As we discussed 
above, the matching procedure is effective to produce new solutions, associated to monodromy data 
for which the connection problem has not yet been studied. Therefore, it is a tool to study the 
classification problem. This classification will be done in another paper. 

A matching procedure, to obtain asymptotic behaviors and monodromy data in the framework 
of the isomonodromy deformation method, was suggested by Its and Novokshcnov in [15], for the 
second and third Painleve equations. The work by Jimbo [17] can be regarded as an implicit 
matching procedure. This method was further developed and used by Kapaev, Kitaev, Andreev, 
and Vardanyan. Here we cite the case of the fifth Painleve equation, in [2] . An analogous matching 
scheme is used in [1], for a different problem (limit PVI — > PV). 

Acknowledgements (May 2006): The author wishes to thank Alexander Kitaev for introducing 
him to the matching procedure and for many discussions. This paper was written in RIMS, Kyoto 
University, supported by the Kyoto Mathematics COE fellowship. 



PART I: Matching Procedure and Results 



2 Matching Procedure 

PVI is the isomonodromy deformation equation of a Fuchsian system of differential equations [18]: 



d5» 

— = A(X,x,6) A(X,x,( 



A (x,9) | A x {x,6) | Ai{x,9) 



A X-x A- 1 



AeC. (4) 



The 2x2 matrices Ai(x, 6) depend on x in such a way that the monodromy of a fundamental solution 
\&(A, x) does not change for small deformations of x. They also depend on the parameters a, /3, 7, S 
of PVI through more elementary parameters 9 — {9q,9 x ,9 1i 9 O0 ) according to the following relations: 

9 1 
-A x := A + Ai+A x = — ^er 3 , 9 X / 0. Eigenvalues (Ai) = ±^ 1 ' i = °' 1,x ' 

a= 1 -(9 oa -ir, -P=\9l 7=1*?, Q-*)=^ (5) 

Here 03 is the Pauli matrix. The equations of monodromy-preserving deformation (Schlesinger 
equations), can be written in Hamiltonian form and reduce to PVI, being the transcendent y(x) 
solution of A(y(x),x, 9)i^2 = 0. Namely: 

V ( X ) = x ( A °) 12 ( 6 ) 

V[ > x p )i2 + (A 1 ) 12 ]-(A 1 ) 12 ' {0) 

The matrices Ai(x, 9), i — 0, x, 1, depend on y(x), dy ^ and j y(x) through rational functions, which 
are given in [18]. In short, we will write Ai = Ai(x). 

The product of the monodromy matrices Mq, M x , Mi of a fundamental matrix solution ^ at 
A = 0, x, 1 respectively, is equal to the monodromy at A = 00. The order of the producs depends on 
the choice of a basis of loops. As a consequence, the following relation must hold: 

cos(7r6» )tr(MiM a; ) + cos(7T0i)tr(M o M a; ) + cos(7r6» :E )tr(MiMo) 

= 2003(7^00) +4008(7^1)008(7^0)008(7^). 
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2.1 Leading Terms of y(x) as a result of Matching 

We present the constructive procedure to obtain the leading terms of a solution y(x), when x — > 0. 
This procedure has been used for the fifth Painleve equation by F.V. Andreev and A.V. Kitaev in 
[2]. An analogous scheme is used in [1], for a different problem. In particular, in [1] the non fuchsian 
singularity in the v]> o;7T -equation (to be introduced below) appears. 

Since we are considering x — > 0, we divide the A-plane into two domains. The "outside" domain 
is defined for A sufficiently big: 

\\\>\x\ SouT , Sout>0. 

Therefore, (4) can be written as: 



(7) 



~dX 



A + A x 
X 



Ax \ - fx 



71=1 



A! 



X 



The "inside" domain is defined for A comparable with x, namely: 

|A|<M 5 '~, S IN >0. 
Therefore, A — > as x — > 0, and we rewrite (4) as: 

A 



~dX 



X X 



n=0 



(8) 



(9) 



(10) 



If the behavior of Aq(x), Ai(x) and A x (x) is sufficiently good, we expect that the higher order 
terms in the series of (8) and (10) are small corrections which can be neglected when x — > 0. If this 
is the case, (8) and (10) reduce respectively to: 
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A T 
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A x 
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where Njn, Nqut are suitable integers. The simplest reduction is to Fuchsian systems: 



d^QUT 

dX 



Aq + A x 
X 



A 1 



A- 1 



^OUT, 



IN 



dX 



Ao A x 
X X — x 



IN- 



(11) 

(12) 

(13) 
(14) 



It is a new feature of this paper that we will use reduced non-fuchsian systems. In the literature, 
the fuchsian reduction has been privileged, but we show that in some relevant cases it cannot be 
used, being the non-fuchsian reduction necessary. 

Generally speaking, we can parameterize the elements of A + A x and Ai of (13) in terms of 6*i, 
the eigenvalues of A n + A x and the eigenvalues #oo of A n + A x + A\. We also need an additional 
unknown function of x. In the same way, we can explicitly parameterize the elements of Aq and A x 
in (14) in terms of 6>o, Xl the eigenvalues of Aq + A x and another additional unknown function of 
x. When the reductions (11) and (12) are non-fuchsian, particular care must be payed. This will 
be explained case by case in the paper. Our purpose is to find the leading term of the unknown 
functions when x — >• 0, in order to determine the critical behavior of Aq(x), A\(x), A x {x) and (6). 

The leading term can be obtained as a result of two facts: 

i) Systems (11) and (12) are isomonodromic. This imposes constraints on the form of the unknown 
functions. Typically, one of them must be constant. 

ii) Two fundamental matrix solutions ^ouriX, x) 7 ^in(X,x) must match in the region of overlap, 
provided this is not empty: 



*OUt{\x) ~ * 7A r(A,x), \x\ So " t < |A| < \x 



(5/ p 







(15) 
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This relation is to be intended in the sense that the leading terms of the local behavior of ^>out 
and vE'/tv for x — > must be equal. This determines a simple relation between the two functions of 
x appearing in A n , A x , A\, A + A x . (15) also implies that Sin < Squt- 



Practically, to fulfill point ii), we will match a fundamental solution of (11) for A 
fundamental solution of (12) when fi := X/x — >• oo, namely with a solution of: 



IN 



dfj, 



IN, 



A 

/i := -. 

x 



0, with a 



(16) 



To summarize, matching two fundamental solutions of the reduced isomonodromic systems (11) 
and (12), we obtain the leading term(s), for x — >• 0, of the entries of the matrices of the original 
system (4). The procedure is algorithmic, no a priori assumption about the behavior being necessary. 

This method is sometimes called coalescence of singularities, because the singularity A = 
and A = x coalesce to produce system (11), while the singularity /i = \ and /i = oo coalesce to 
produce system (16). Coalescence of singularities was first used by M. Jimbo in [17] to compute 
the monodromy matrices of (4) for a class of solutions of (PVI) with leading term y(x) ~ a x 1 ~ a , 
< Kcr < 1. 



2.2 Computation of the Monodromy Data 

Let W be a fundamental matrix solution of (4), and let M , M x , M l7 M^ be its monodromy matrices 
at A = 0, x, 1, oo respectively (M^ is the product of M , M X1 M l7 the order depending on the choice 
of a basis of loops). As a consequence of isomonodromicity, there exists a fundamental solution 
^out of (11) such that 



M? UT = Mi , M° UT = M, 



OO ) 



where Mi UT and M^ UT are the monodromy matrices of ^out at A = l,oo. Moreover, Mq UT — 
MqM x or M x Mq, depending on the order of loops. A detailed proof of these facts can be found in 
[9]. There also exists a fundamental solution * JAr of (12) such that: 

M /w = M , Mi N = M x , 

where Mq N and M X N are the monodromy matrices of ^in at A = 0,x. 

The method of coalescence of singularities is useful when the monodromy of the reduced systems 
(11), (12) can be explicitly computed. This is the case when the reduction is fuchsian (namely 
(13), (14)), because fuchsian systems with three singular points are equivalent to a Gauss hyper- 
geometric equation (see Appendix 1). For the non- fuchsian reduction, in general we can compute 
the monodromy when (11), (12) are solvable in terms of special or elementary functions. This will 
be discussed case by case in the paper. 

In order for this procedure to work, not only ^out and ^in must match with each other, as in 
subsection 2.1, but also ^out must match with a fundamental matrix solution <J> of (4) in a domain 
of the A plane, and <J in must match with the same in another domain of the A plane. 

The standard choice of <J> is as follows: 

[ / + (i)] A-^A*-, A^oo; 

Mx)[l + 0(X)] X^X R °C Q , A^O; 

<p x (x)[l + 0(X-x)] {X-x)^{X-x) R *C x , X^x; 

{ Vi(z)[/ + 0(A-1)] (A-l)^(A-l) fll Ci, A-H; 

Here ipo(x), ^Px(x), i>i(x) are the diagonalizing matrices of A (x), A\(x), A x (x) respectively. They 
are defined by multiplication to the right by arbitrary diagonal matrices, possibly depending on x. 



tf (A) = < 



(17) 
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C K , k — oo, 0, x, 1, are invertible connection matrices, independent of x [18]. Each R Kl k = oo, 0, x, 1, 
is also independent of x, and: 



r k = o if e K £ Z, R K = I 



* 





* 



, if 6 K > integer 
, if 6 K < integer 



If 9i = 0, i = 0, x, 1, then i?j is to be considered the Jordan form [J? ^ ) of A.j . If 9^ = 0, Rco = 0. 



x 0, 

Note that for the loop A i-> Ae 27 ™, |A| > max{l, x|}, we immediately compute the monodromy at 
infinity: 

Mqo = cxp{— inOoo} exp{27rii?oo}. 

Let ^out and ^jjv be the solutions of (11) and (12) matching as in (15). We explain how they 
are matched with (17). 

(*) Matching * o ^out- 

A = oo is a fuchsian singularity of (11), with residue —Aoo/X. Therefore, we can always find a 
fundamental matrix solution with behavior: 



.j t Match 
^OUT 



A-^A^, A^ 



oo. 



This solution matches with <J/. Also A = 1 is a fuchsian singularity of (11). Therefore, we have: 
^atch = ^,out {x ^ i + 0{x _ 1 ^ (A-i)^(A-l)^Cf t7T , A-H; 

Here Cf" 77 is a suitable connection matrix. ■0? C/T ( a; ) is the matrix that diagonalizes the leading 
terms of A\{x). Therefore, ipi(x) ~ \pi UT (x) for x — > 0. As a consequence of isomonodromicity, i?i 
is the same of 

As a consequence of the matching ^ <H> ^ M ut'\ the monodromy of ^ at A = 1 is: 
Mi = Cr 1 exp{i7r6»icr 3 } exp{27rii?i}Ci, with Ci = C? UT . 

We finally need an invertible connection matrix Cout to connect ^ M ff^ h with the solution ^out 
appearing in (15). Namely, ^ M {%l h = ^ olJT c OUT . 

(*) Matching * o */at: 

As a consequence of the matching VP -B- ^out* 1 ' we have to choose the IN-solution which 
matches with V M ffi h . This is ^ Match . = ^, inCo[jt 

Now, A = 0,x are fuchsian singularities of (12). Therefore: 

^ N (x)[l + 0(\)] \^^\ R «C£ N , A^0; 

^ X N {x)[l + 0{\ - x)} {\-x)^{\-x) R *Cl N , A^x; 

The above hold for fixed small x ^ 0. Here C^ N and C£ N are suitable connection matrices. ipQ N (x) 
and il) x {x) IN are diagonalizing matrices of the leading terms of Aq(x) and A x (x). For x — > they 
match with tpo(x) and ^(z) of W in (17). On the other hand, as a consequence of isomonodromicity, 
the matrices i?o and R x are the same of 

By virtue of the matching ^ <H- \E , ^ tc/l , the connection matrices C and C x coincide with the 
x- independent connection matrices Cq N , C 7Ar respectively. As a result, we obtain the monodromy 
matrices for <]/: 

M - Co' 1 eMi^oos} exp{2 7 ri J Ro}Co, C = C 7Ar , 
M x = C x - X exp{in0 x a 3 } cxp{2niR x }C x , C x = C ! X N . 

Our reduction is useful if the connection matrices Cf UT , Cq N , C 7JV can be computed explicitly. 
This is possible for the fuchsian reduced systems (13), (14). For non-fuchsian reduced systems, we 
will discuss the computability case by case. 



vr,Match 

W IN — 
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3 Results 

In the following, it is understood that i^O inside a sector. Namely, arg(x) is bounded. 
3.1 Critical Behaviors: Result R3). 

The novelty of this paper is that the matching procedure is applied to non-fuchsian systems (11) 
and (12). As a result, we obtain all the solutions that admit a Taylor expansion 

oo 

y(x) = b Q + bix + b 2 x 2 + ... = ^ b n x n , x -+ 0. 

Precisely, we obtain the representative solutions of three equivalence classes, the equivalence relation 
being the birational transformations [24] of Appendix 3 and formula (22). Our result is the following. 

Theorem 1 The solutions of (PVI) with Taylor expansion at x = are divided into four equivalent 
classes (one being that of singular solutions y = 0,1, x). The representatives can be chosen as 
follows: 

1) Singular solution y = 1. 

2) 0oo ^ 1, 0i - 0oo <£ Z [representative of 1 ±0 OC ^ Z]: 

/ \ 0\ - 0oc + 1 . 0\[{0\ - 0oo){0\ - 0oo + 2) + 9 2 - 6q] ^ ulna a a \ n n o\ 

y{x) = ^iLr + 2(^-1)^-^(^-^-2) x 5 00 x ' (18) 

The coefficients are certain rational functions of 0q, 0oc, 0o, 0x- 

3) 0i=0oo^ 1, 6*o = ±0x [representative of 1 ±6 oo e Z, x ±0 Q e Z]: 

y(x) = h ax + \ 6 n (a; 6> , (19) 

T/ie coefficients are certain rational functions of 0q, 0oo and a parameter a 6 C. 

4) #oo = 1, 0i = [representative of 6>i ± 0^ e Z, fl^ e Z\{0}]: 

_ 00 
y{x)=a+^(l + 01-01) x +J2b n (a;0 ;6 x )x n . (20) 

n=2 

The coefficients are certain rational functions of n , X and a parameter a E C. 

The monodromy data associated to the above solutions is given in theorem 3. 
The symmetry 0\ n> —0\, which leaves (PVI) invariant, transforms (18) into: 

y(x) ~ x -i + 2(1-0^ + 8^ + 0,-2) *+2J>»i-Woo,e ,e x )z . (2i) 

Here ^ ^ 1, 9\ + #oo Z. The coefficients 6„ are the same of (18). 

The convergence of the Taylor series can be proved by a Briot-Bouquet like argument. This will 
not be done here, for reasons of space. The reader can find the general procedure in [14] and an 
application to the fifth Painleve equation in [20] 

Comments: 

1) Characterization of solutions y(x) = Y^=ob n x n , bo ^ 0. 

(a) There always exists one solution (18) when 0\ — X $ Z; there always exists one solution 
(21) when 0i + 0oc ^ Z. The coefficients b n depend rationally on Kl k = 0, x, 1, oo. (b) There is a 
one-parameter family of solutions equivalent to (19), when 0i ±0^ £ Z and 0q ± 0x has a particular 
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integer value. The coefficients b n depend rationally on a complex parameter a and 6*00,6*0. ( c ) 
Finally, there is a one-parameter family of solutions equivalent to (20), when 9\ ± #00 G Z, and 6*00 
has a particular integer value; the coefficients b n depend rationally on a complex parameter a and 
8 ,9 X . The singular solutions y = 0,l,x are possibly obtained by birational transformations of (18), 
(19), (20). 

The coefficients b n can always be computed recursively by direct substitution into (PVI). We 
will clarify these facts by some examples in Appendix 4. 

2) Characterization of solutions y(x) = Y^=i^nX n , 61 ^ 0. 

These solutions are obtained from those of theorem 1 by the symmetry. 

9 X ^6 U o i->0oo-l, 0i^6 x , flooh-^o + l; y{x) H- (22) 

2/0*0 

The solutions obtained from the singular solution y = 1 and (18), (19), (20) are respectively: 

1) Singular solution y(x) = x. 

2) O ^0,6 ±0 X £ Z: 



3) O + 0* = 1, O 0, 0i = ±(0co - 1): 

CO 

y(a;) = 0oa; + a x 2 + ^ 6„(a; O , 0oo)a; n . (24) 

n— 3 



4) 9 X = O = 0. 

y^) = ax + ^-^(0? - (0oc - l) 2 - l)x 2 + ^ b n (a; 9 1 ,9 0O )x n . (25) 

n— 3 



(a) (PVI) has always one or both solutions (23) when 9$±9 X £ Z. Also when O + 9 X (or O — X ) 
is integer, (PVI) has a solution (23) corresponding to 0o — 9 X not integer (or 0o + 9 X not integer), 
(b) When O + 9 X or O — 9 X is integer, (PVI) has a 1-parameter family of solutions equivalent (by 
birational transformations) to (24); this family exists provided that 9\ ± 0oo has a particular integer 
value, (c) When O + 9 X or O — 9 X is integer and O has a particular integer value, there is a one 
parameter family of solutions equivalent to (25). 

3) (PVI) has a one-parameter family of solutions of the type: 

00 

y(x) = y {x) + (x) ax" + y 2 (x) (ax") 2 + ... = ^ y N (x) (ax") N , x -> 0; (26) 

where the parameter is a £ C, and the j/Af(x)'s are Taylor series: 

00 

Un(x) = ^^^(01,000,00,02;) X k , X^O. 

Either yo(x) is (21) and uj = ±(01+0^ — 1), or yo(x) is (18) and w = ±(0oo — 0i — 1) • The conditions 
|SRw| < 1, w / hold. The coefficients &fc,Ar(0i, 0oo, 0o, 0^) ar e certain rational functions that can 
be recursively determined by direct substitution into (PVI). These solutions are the immages of 
solutions (28) and (29) respectively, through the symmetry (22). Solutions (28) and (29) are a 
sub-case of theorem 2, obtained by the matching procedure. 

Taylor solutions (18), (21) are a special case of (26), when the parameter is zero. Solutions (19) 
and (20) - and their images by symmetry - are one parameters families of type (26), in non generic 
cases when w 6 Z. 
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Further study of one-parameter solutions, including non-generic cases when 9 V and/or some sum 
of two #„'s are integer (including logarithmic one parameter families), will be presented in another 
paper devoted to the general classification problem. 

4) Solutions (18) and the equivalent solutions (21), (23) were also derived in [19] by substitution 
of a Taylor expansion in (PVI). The corresponding monodromy was computed by coalescence of 
singularities of a Heun's type (scalar) equation. 



3.2 Critical Behaviors: Results Rl), R2). 

We now consider cases when (4) can be reduced to the fuchsian systems (13) and (14). Let a be a 
complex number defined, up to sign, by: 

tr {M M X ) = 2 cosher), \$ta\ < 1. 

Actually, ±<r/2 are the eigenvalues of lim^o^o + A X ). The matching procedure yields the following 
result. 

Theorem 2 Let r E C and a be as above, with the restriction \$ta\ < 1. (PVI) has a family of 
solutions depending on the two parameters r, a. The leading terms of the critical behavior for x — > 
may be parametrized as follows: 



For a ^ 0: 



( l [v 2 -(e„+e x ) 2 ][(e -e x ) 2 -* 2 ] 



y(x) ~ < 



16<T 3 



x 1+<T , 



x ^iA sin(ioTna; 
In the above formulae, r ^ and 



"ST* - 



■}• 



if Kcr > 0; 
if 5Rcr < 0; 
if SRcr = 0. 



(27) 



i In — -, A := 
a A 



J 2 I 2a 2 



For special values of a j^O: 



For a = 0: 



y{x) 



y{x) 



x =F 



0q + Ox 0q 



0o 



0q - o 3 



x =F 



0o — X 



,1+cr 



y(x) ~ < 



ei-ei 



In x + ^ffrfsi 



4r+26 
01-01 



a = ±(0 O + 9 X ) + 0, 
a = ±(0„ - 9 X ) + 0. 



S— 



x (r ± 9 Inx), 



0o i= 



(28) 
(29) 

(30) 



Comments: 

1) r can be computed as a function of the monodromy data. See (37) and comments there. The 
sign (branch) of the two quare roots appearing in <f) and A is the same, x — > in a sector of width 
less then 2ir. 

2) Sub-cases of theorem 2. 

i) When a ^ 0, the result of the theorem includes the sub-cases (28) and (29). If r = 0, 9 n ^ 0, 
#o ± X £ Z, direct substitution into (PVI) gives the two Taylor expansions (23). 
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If r ^ 0, (28) and (29) are a 1-parameter family, with the restriction |3?cr| < 1. The symmetry 
(22) transforms them into the solutions (26), the leading terms being respectively: 



y{x) 



y(x) 



&OD +6*1 — 1 
0oo — 1 



1 ± 



#oo — 1 



1 ± 



- 1 



with the restriction |SRw| < 1 . 

ii) The case a = includes the sub-case y(x) 
substitution in (PVI) we obtain a series: 



W = ±(000 + 01-1)^0, 
U = ±(0oo - 01 - 1) ^ 0, 

rx, which occurs for O = X , O = 0. By direct 



y{x) 



r x 



E 

n=3 



b n (r,0 1 ,0 oo )x n , O =0 X = O, r^0,l. 



This is again solution (25). Note that for r = 0, 1 we have the singular solutions y = 0, y = 1. Note 
also that the special sub-sub-case 8o = X = Qi = has applications in the theory of semi-simple 
Frobenius manifolds of dimension three [7] [10]. 

3) The first two behaviors (27) were derived in [17]. The existence of such solutions was proved 
by assuming that the matrices Aq, A x , A\ have a certain critical behavior for x — > 0, and proving 
that such matrices solve the Schlesinger equations. Then, the monodromy data were computed by 
coalescence of singularity. These solutions where further studied in [8], [9], [11], [5]. We show that 
these solutions can be obtained without any assumption by the matching procedure, together with 
the solutions (30) and the first of (27), which do not appear in [17]. 

4) The class of solutions (27) was enlarged in [26] and [11], to the values a G C, a £ (-co, 0]U[1, +oo). 
When Sftcr > 1 or 5ftcr < 0, the critical behavior is like the first of (27), and it holds for x — > in 
a spiral-shaped domain in the universal covering of a punctured neighborhood of a; = 0, along a 
paths joining a point xo to x = 0. Along special paths which approach the movable poles, these 
solution may have behavior y(x) ~ sin - \nx+ip(x, r)) , where if(x, r) is a phase depending on the 
parameter r. The transformation a ±a + 27V, N e Z, leaves the identity ti(M M x ) = 2cos(7rcr) 
invariant. Its effect on the solutions is studied in [11]. As a result, one can reduce to the values 
< Sftcr < 1, g 7^ 0, 1. We cannot enter into more details here. The reader may find a synthetic 
description of these results in the review paper [12]. 

5) Solutions with expansion: 



y(x) = x(A x + Si Inx + Ci In 2 x + Di In 3 x + ...) + x 2 (A 2 + B 2 hix + ...) + a 
are all included it theorems 1 and 2. Actually, only the following cases are possible: 

"° -x + 0(x 2 ) [Taylor expansion], 



0. 



y(x) = < 



x (^J+Bilna; + 
x (A 1 ±6 \nx)+x 2 (...) + 



and 6n — ±9 X 



(31) 



Ai and B\ are parameters. We see that the higher orders in (30) are 0(x 2 ln m x), for some integer 
m > 0. 



6) The symmetry (22) applied to solutions (30) gives: 

4r + 2(6 



namely: 



and 



y(x) ~ 4 [(9j - {6 X - l) 2 ) 



y{x) 



Ins — 



-1) 



[6\ - (0oo - l) 2 ] ln 2 a; 



- (^oc - l) 2 

8r + 4(000 - 



4(0 O 



1) 



91 - (0oo - l) s 



1 

In a; 



O 



\n 2 x 



(32) 



y{x) = 



±i 



(floe - 1) lnx 



1 + 



(0oo-l) lnx 



+ 



ln 2 .x 



0oo +01 — 1- 
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The higer orders 0(1/ In 2 x) include powers x n (lnx) ±m . The so called Chazy solutions, studied in 
[21] for the special case 9 n = 9 X = 9\ = 0, 9^ = —1, have the behavior (32). 

7) When this paper was completed, I received a communication by the first author of [6] . In [6] it is 
proved that (PVI) has solutions with expansion at x — oo, or x = 0, of the form y — c r x r + ~^2 S c s x s , 
c r G C. The c s 's are either complex constants or polynomials in lnx. r and s are integer or complex. 
If r is complex, the restriction 5Rr e (0, 1) holds. The method used in [6] is a power geometry 
technique. The connection problem and the characterization of the associated monodromy data are 
not studied. 

8) When this paper was already accepted for publication, I received a private communication about 
a recent work [25], on the asymptotics of the real solutions of (PVI). The asymptotic behaviors 
obtained in [25] are of the type of our theorem 2, namely (27) and the first behavior in (30). The 
tool used is a method of successive approximations. So, the results are local, and the connection 
problem is not studied. Moreover, some genericity conditions on the coefficients of PVI seem to be 
necessary (so, for example, the second solution in (30) cannot be obtained). 

3.3 Monodromy: Result R4). 

In this paper, we computed the monodromy for the Taylor-expanded solutions, which correspond to 
non-fuchsian reductions of system (4). Because of the symmetries of (PVI), we can limit ourselves 
to the monodromy data for the representative solutions (18), (19) and (20). 

Theorem 3 a) Let 9 K ^ Z, n — 0,1, x, oo. A representation for the monodromy matrices of the 
solution (18) is: 

M = C 0oo e?q>{iir9 a o- 3 } C^, 
M x = Cooo Cqi cxp{iir9 x cr 3 } C i C^. 
Mi = exp{-i7r#icr 3 }, M M = cxp{-iTr9 00 a 3 }. 
The matrices Cqoo and Cqi are: 



- r (i+^-^)r(i+e )e I f f>a+^+»^-^ r(i+^-%L)r(i-0o)e < i [e *- e ° +e °°- 911 
r(^ + %^-%+i)r(%-^+2i-^ + i) r(-%-^-% + ^+i)r(^-% + ^-%+i) 
r(%- -■^-l)^(l+0o)e i '^' [^,o+£)x+£)l ~ ^,ool p^Qoo »\ i)r(i e.-jc^^'""^ 1 ~ e °° ] 

L L \ 2 ~*~ 2 ^ 2 2 J L \ 2 2" r 2 2 J L \ 2 2 2 ^ 2 ) L \ 2 2 ' 2 2 ) -I 

(33) 



Coi 



r(-9 x )r(i+e ) r(-e x )r(i-e ) 

r(^-^+^-^+i)r(^-Y+Y-^) r(-^-^_2 f , + «j.; i )r(-^-^-^+Y) 

r(6) x )r(i+e ) r(e x )r(i-e ) 

r(%+^+%-^)r(^ + ^+^-%+i) r(^-$ + %^)r(%L-£+£-^+i) 



(34) 



T/ie subgroup generated by M M X and Mi is reducible. As for the solution (21), we just need to 
change 9 X i-> —8\. 

b) It is convenient to re-parameterize the solution (19) by introducing a parameter s through the 
equality: 

= ^(25 + ^ + 1) 

a 2(0^ - 1) ' 
Let 9 X , ^ Z. Then, a representation for the monodromy group is: 

M = G exp{iir9 x a 3 } Mi = exp{ —^9^ a 3 ] 

M x = G eTip{-iir9 x a 3 } G~ x , M^ = cxp{-i7r6»ooCT3} 
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In particular, Mi — M^, M$M X = I. We can choose G as follows: 

r-f 1 1 

\ r r 

Conversely, we may express s as a funcition of the monodromy data: 

O x [2 cos(7r(0oo + 9 X )) - tr(Mi M )] 



s = 



2 [003(^(6*00 - e x )) - cosinie^ + e x ))] 



c) We re-parameterize solution (20) introducing a new parameter s defined by a =: (1 — s) 1 . 
Let 9 , 9 X £ Z. Then, a monodromy representation for the solutions (20) is: 



M = (Ccoo) 1 cxp{z7r6» (73} C^o, 



M n 



-1 
2m (1 - s) -1 



M x = (Cooo) 1 (C i) 1 exp{iir6 x (j 3 } CmCooo, M x = 



1 

2nis 1 



where and Cq\ are (35) and (36) given below. Conversely, we may express s as a function of 
the monodromy data: 

tr(MiMo) - 2cos(ir0 o ) (C , oo0 ) 2 i 



s = 



47T sin(7r6'o) 



T/ie matrices C^q and Cqi are: 



ooO J22 



CooO 



r(-»-^_i)r(-^+^- 



r(i-e )e" 



■ r»0 fx 3i 

1 9 9 1 ; 



r(-8 ) £ 



,r«0 ,«I .8, 



Jo _<> 



r(flo) 



(35) 



r(-6) x )r(i+e ) 



Ci 



01 



r(-e x )r(i-e ) 



r(8 3; )r(i-6io) 



The conditions 9 K ^ Z can be eliminated, and the computations can be repeated without con- 
ceptual changes, but with different results. 

In the above theorem, the subgroups generated by M M X and Mi are reducible. This char- 
acterizes the monodromy associated to solutions which have a Taylor series at x = 0. The same 
characterization at x = 1 involves the subgroup generated by M\M X and Mq. At x = oo, it in- 
volves the subgroup generated by M Mi and M x . 2 In another paper, we will consider again this 
characterization, together with the general problem of classification. 

Let us define again a by tr(MoM x ) = 2 cos7rcr. Then, in case a), a — ±{9 1 —9 O0 ) [and ±(#i +#00) 
for the change 9 X -9{\. In case b), tr(M M x ) = 2 and a = 0. In case c), tr(M M x ) = -2, a = ±1. 
The matching procedure is effective to produce solutions corresponding to monodromy data for which 
the connection problem is so far not well studied, such as the case tr(MjMj) = —2 (see [12]). 3 

Note: Also the 1-parameter solutions (28) (29) and the second solution in (30) are characterized by 
a reducible subgroup generated by M , M x . 

2 In the appendix of [11], the reader may find explanations about how to obtain results at x = l,oo from the 
results at x = 

3 Here I remark that the formula (1.30), page 1293, of my paper [11] is wrong. The correct one is tr(M;Mj) 
(-00, -2]. In [11] the connection problem is solved for tr(MiMj) ^ ±2. The case tr(MjMj) = 2 yields (30). For the 
special choice of the parameters 80 = 8 X = d\ = 0, it was studied in [8] and [9] (no logarithmic terms appear in such 
a special case). The result (30) for the general (PVI), corresponding to tv(MoM x ) = 2, appears in the present paper 
for the first time. 
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Comments. 



1) The monodromy group for the solutions (23) was derived also in [19], by confluence of singularities 
of scalar equations (including a Heun's type equation). The result is equivalent to that in point a) 
of the above theorem. 

2) The computation of the monodromy group of the fuchsian systems (13) and (14) is quite clear 
[17] [8] [11] [5]. It allows to express the parameter r of (27), (28), (29) and (30) as a function of the 
monodromy data. Though the computation for (28), (29) and (30) does not appear in the literature, 
the procedure is clear (see section 4.8), so we do not repeat it. We just report the result for (27), 
which can be found in [17] [11] [5]: 



where 



and: 



{e -e x + cj){e + e x -a){e oo + e 1 



_£) I 
F' 



F := 



U :-- 



4(7(6*00 + 6>i + cr) 

r(i + a) 2 r (i(flo + e x - a) + 1) r (\{e x -e a -a) + i) 

r(l - afT (±(0 O + 6 X + a) + l) T (\{6 X - O + a) + l) 

r (±(0oo + 0i - a) + 1) r (±(0i - 0QQ - a) + 1) v 
X r (±(000 + 0! + a) + 1) r (±(0i - 0oo + a) + 1) u ' 

■ sin(7T(j)tr(MiM a ;) — cos(tt8 x ) 003(^000) — cos(7T0 o ) cos(7T0i) 



+ 



(37) 



+- sin(7ra)tr(MoMi) + cos(7r6 x ) cos(7r#i) + cos(7T0oo) cos(7r# ) 

7T 7T 7T 7T 

V := 4sin-(0 o + 8 X - a) sin -(0 O - 0^ + cr) sin-(0oo + 0i - a) sin - (0^ 



+ a). 



The above formula was computed with the assumption that cr±(0o + 0a;), cr ± (0o — 6 X ), cr±(0i + 0OO ) , 
cr ± (0i — 0oo) are not even integers. 4 

3) Reducible Monodromy. The monodromy groups in theorem 3 are not reducible, but they 
have a reducible subgroup. If the entire group itself is completely reducible, all the Painlevee 
transcendents are known. Solutions of (PVI) corresponding to a reducible monodromy were found 
in [13]). We summarize the results: 

Proposition 1 All the solutions of (PVI) corresponding to a reducible monodromy group are equiv- 
alent by birational canonical transformations to the following one-parameter family of solutions, with 
0oo + 0i + 0o + &x = 0: 



h + Goo — 1 + x(l + 9 X ) 1 x (1 — x) du(x;a) 



- 1 



0oo — 1 u(x; a) 



dx 



(38) 



where u(x;a) = U\{x) + au2(x); a £ C, U\(x) and U2{x) are linear independent solutions of the 
hyper- geometric equation: 



ffi~ It fill 

< l - x )^2 + it 2 - ( d °° + - ( 4 - 6 °° + W ^ 



(2 - oc )(i + e x )u = o 



The monodromy matrices are 



in 



^=(o M -=U -th *=U -* 

The parameter a does not appear in the monodromy. 

4 In [11] there is a missprint in formula (A. 30), which must be re-calculated. In [17], in formula (1.8) at the bottom 
of page 1141, the last sign is ±<r instead of =f<t. 
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Remark: The rational solutions of (PVI) are a special case of the above proposition. They were 
studied in [22]. Up to canonical birational transformations, they are realized for + 6\ + + 6 X = 
and: 

a , / \ Oqq + 81 x - 1 

y[x> 0co x(i + e 1 )-(e 1 + e oc )> 



-2 : 



y{x) 



(2 



h x) -2 



1\ x 



(i-9 0O )(2-{e 0O + e 1 ) + e 1 x ) 



The computation of the expansion at x — of (38) is just a consequence of the expansions of ui (x) 
and v,2(x). The reader can find by himself a behavior y ~ x(r(a)±6 x ln(x)) for di+d^ = 9 + 8 x = 0, 
namely a sub-case of the second solution in (30). For 81+800 Z, we find behaviors of the type 
(26) (and (18), (23) for a = 0). 



PART II — Derivation Results: Fuchsian Reduction 



4 Fuchsian Case 

Let x — > 0. The reduction to the fuchsian systems (13) is possible if in the domain (7) we have: 



|(A) + A x )ij\ > (A x )ij - , namely: \{A + A x ) tJ \ > \(A x ) l3 x 1 ~ s ° UT \ . (39) 
Let us denote with Ai the leading term of the matrix Ai, i — 0,x, 1. We can substitute (13) with: 

d^QUT 



dX 



A a + A x A 1 



A 



A- 1 



^OUT 



(40) 



We suppose that 8^ ^ 0. This is not a loss in generality, because 800 = is equivalent to 8^ = 2. 

Lemma 1 If the approximation (13) is possible, then A + A x has eigenvalues ±£ G C independent 
of x, defined (up to sign and addition of an integer) by tr(M x M ) = 2cos(7rer). Let n 6 C, ri ^ 0. 
For 8oo ^ 0, the leading terms are: 



V 160^ ri 



-n 

tT 2_g2 ; _g2 I , 
40 ro / 



and 



A +A x 



k 2 -(8i-8oc) 2 ]k 2 -(ei+8oc) 2 ] j_ 



(41) 



(42) 



Proof: Observe that ti(A + A x ) = tr(A + A x ) = 0, thus, for any x, A + A x has eigenvalues 
of opposite sign, that we denote ±cr(x)/2. Then, we recall that x is a monodromy preserving 
deformation, therefore the monodromy matrices of (40) are independent of x. At A = 0, l,oo they 
are: 



AlOUT _ I M X M A/ rOUT _ M AjOUT _ , , 

Thus, det(Mo UT ) — 1, because dct(Af 2; )=dct(Mo) = 1. Therefore, there exists a constant matrix 
D and a complex constant number a such that: 

diag(exp{— Z7r cr}, exp{z7rtr}), 

or ( # ±1 I , ,€Z 
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We conclude that a(x) = a. We also have tT(M DUT ) = 2 cosher). 
Now consider the gauge: 



$i:=A 2 (A-1) ^ * oc/T . ___ 

(LA 



A + A x - f ^ - f 



A 



A- 1 



$1 



(43) 



We can identify A + A x — ^ and Ai ~ ^- with _B and -Bi of Proposition 2 in Appendix 1, case 

(98) , with a=^f + ^ + f,fc = + \ + § , c = a. □ 

Remark: n may be a function of x. If the monodromy of system (40) depends on n, then n is a 
constant independent of x. This is the case here, but we do not prove it for reasons of space. Sec 
the references in Part I. 

Lemma 1 (and Lemma 3 which follows) includes all cases (98) — (102) for system (43). Cases 

(99) — (102) are obtained substituting a = — (#oc+#i), #oo — #i> #oc + #i, #1— #00 respectively. For all 
the computations which follow, involving system (40) or (43), we note that the hypothesis S M / 
excludes cases (103), (104) and the Jordan cases (105)-(107). 

The reduction to the fuchsian system (14) is possible for x — > in the domain (9) if: 



(Ao)ij , (A x )i 



> |(Ai)y | , namely: 



(Aa+A^ij 



A A — x 

We can rewrite (14) using just the leading terms of the matrices: 



» \(Ai)ij\ 



IN 



dX 



A 



A A — x 



IN, 



Then, we re-scale A and consider the following system: 

d^IN [Ao, A* 



d/j, \ /j, /j, — 1 I 
We know that there exists a matrix K (x) such that: 



IN, 



A 

/1 := - 









I— 1 


(0 -§) 


or 


s 





(44) 



(45) 



Let Ai := Kq AiKo, i = 0,x. By a gauge transformation, we get the system: 



* /A r =: K (x) * , 
Lemma 2 Let reC,r/0. If a ^ 0, we have: 

A Q = 



d*o 
dfi 



An 



H (J.-1 



l-el+a 2 



4a 

g 2 -(e a -e :E ) 2 }[a 2 -{e +e :L ) 2 } j_ 



ro 



4a 



4cr 



k 2 -(eo-e 3; ) 2 ]k 2 -(eo+9 3; ) 2 ] 1 cr 2 +sj-e 2 a 



(46) 

(47) 
(48) 



16<T 2 



4(7 



Proof: We do a gauge transformation: 



$0 :=A* 1) *o, -j— 



d£o 



A - 



+ 



(49) 
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We identify A — f , A x — f with B and Bi in the Appendix 1, Proposition 2, case (98), with 
a=% + %- + % + §,c = 6> . □ 

Remark: If the monodromy of the system (46) depends on r, then r is a constant independent of x. 
This is the case here. 

The Lemma 2 (and lemma 4 which follows) includes also the cases (99)-(102) for the system 
(49). These cases correspond respectively to the values a = 9 + X , —9 — 9 X , 9 X — 6o, 9 — X , with 

#0 / ±9 X . 

4.1 Matching for a £ Z and proof of (27) 

We match "if out and ^jjv in the intersection of the "outside" and "inside" domains, namely the 
region Ixl 50177, < |A| < |a;| 5/N , x — V 0. As a consequence, we obtain the leading term of y(x). 

Lemma 3 If a £ Z and 9^ ^ 0, system (40) has a fundamental matrix solution if? out W with the 
following behavior at A = 0: 

^ouT = f^G n \ n ( Aa A), Go = f (e^-lf-el ) ■ 

n=0 V / V 4eoon 40 TO n / 

G„ are matrices which depend rationally on 9^, 9\, a, r\. The series is convergent for |A| < 1. 

Proof: It is an immediate consequence of the standard theory of linear systems of fuchsian differential 
equations. □ 

Lemma 4 If a $ Z, system (46) has a fundamental matrix solution with the following behavior at 
fi = oo: 

jJL* 



*o(/i) 



n=l 



jU" 



where I is the identity matrix, K n are matrices which depend rationally on 9 , 9 X , a, r. The series 
is convergent for \fj,\ > 1. 

Proof: It is a consequence of the standard theory of systems of fuchsian equations. □ 
The matching relation *i(A) ~ K (x)$? (A/a;), \x\ s ° UT < |A| < \x\ SlN , x ->• 0, is: 

This gives the result: 



A"t ) V A"* ) \ 



1 1 \ I 12 

40 oo r 1 46 00 r 1 



a;-2 



We compute the matrices A (x) = Ko(x)A K (x) , A^x) = ifo(£)^-^o(x) 1 making use of 
Lemma 2. We obtain: 

= G I ((T+ g x _g o)(<J+ g x+ g o)(CT ' T _g x+ g o)(CT _g x _ t)o) eg-e^+T 2 J G ° 

\ 16^ X 4cr / 

A x {x) = Go ( (CT+ e x _9 )( CT+9x+ e o)(CT % x+ e o)(CT _e x _9 o) CT 2 +e 2_ 9 2 ) G ° 

V 16<r 2 r X 4cr / 

This result shows that the matrix elements of Aq and A x diverge as |a;| _ l 3?CT l when x — > inside a 
sector (i.e. for |arg(x)| bounded). In particular, we find (Ai)i2 = —r\ and 

(Aj)l2 = 7" 16^ X + ri 2^ V X ■ 
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The above are enough to compute the leading term(s) of y(x) from the formula: 



, , x(A ) 12 x(A ) 12 

V[X) ~ x[(A ) 12 + (^) 12 ] - (A!) 12 - (A,) 12 

Thus: 



(Ai)i 2 



! -1 



, . (A )l2 

y(x) ~ —x 



1 [a 2 -(e„ + ^) 2 ][(e -^) 2 -a 2 ] !_ CT 0,i-0 2 +a 2 r 1+CT 

" ic3 * H x x 

r lbcr A la 1 a 



(50) 



(51) 



(Ai)i 2 

We have ignored 1/[1 — a;(l + (^4.0)12/(^.1)12)] because condition (39) is equivalent to: |a; 1_ ' 5l± ' T | — > 
for x -> 0, which implies that x (^0)12/(^1)12 ~ \x 1±a \ -> 0. Therefore, 1/[1 - x(l + 

(io)i 2 /(ii)i 2 )] = 1/(1 + 0(af)) = 1 + O(x). 

If Sftcr ^ 0, the leading term of (51) is certainly correct, but some higher order corrections may 
be bigger than the next two terms of (51). If 3?cr = 0, the three terms of (51) are of the same order, 
and their combination gives the trigonometric expression in theorem 2. 

4.2 Range of a 

Conditions (39) and (44) must be verified. Let C denote a non zero constant. We suppose that 
x — V inside a sector with center on x = 0. Then: 

Condition (44) is \x\~ SlN > C S IN > 0. 

Condition (39) is C> \ x \-\*"\+i-Sout ^ j^j < x _ £ ot/T . 
The last condition implies that |5Rcr| < 1. We also conclude that < Sin < Sour < 1- 

4.3 Leading term for a = ±(0„ + 9 X ), ±(0 O - 0*) 7^ 0. Proof of (28) and (29) 

Formula (51) holds for any a ^ such that |5ftcr < 1. However, we cannot naively substitute the 
value of a — ±(6> + X ), ± (#o — X ), f° r which the coefficient of x 1 ~ a vanishes. This is because only 
the leading term is certainly correct, and it may be the term in x 1 ~ a . Therefore, here we briefly 
give the explicit derivation of (28) and (29), using cases (99)-(102) for system (49). 

Case (99), a = 0: This is the case 7 = ^ + ^/ 0. The matching procedure does not change. From 
(99) we compute: 

( AO r ™<3 \ / Ox r ^ \ 

40 Go " Ax=Ga (0 -4 J Go " 

This implies that (A a )i 2 = r\ (yj^fg 8o + g x a ^j , while (Ai)i 2 = —r\ as in the generic case. 

Therefore, 

y(x) ~ — x — — x T . 

0q + X 8q + X 

It is interesting to note that for Sfter > we have y(x) ~ 6 /(8 + X ) x. Such a behavior is what one 
would naively expect from the generic behavior (51) when a = 0. 

Case (100), b = 0, is a = -6 - 6 X ^ 0. Case (101), a = c, is a = 6 X - O . Case (102), b = c, is 
(j = $0 — ^2: • Proceeding as above, we find (28) and (29). 

Remark: If we substitute y = b\x + b 2 x 2 + b^x 3 + ... into (PVI) we find all the coefficients b n by 
identifying equal powers of x. The result is (23). We need to assume that #0 ± X is not integer or 
zero. 
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4.4 Matching for a = 0. Proof of (30) 

4.5 Case O ± 6 X ^ 

Lemma 5 Let ri e C, ri ^ 0. TTie matrices of system (40) are: 

/ floo 2 + 8l 



^1 



40^" "ri \ - r 

?_fl 2 l 2 n 2 , a 2 \ , Aq+ A x 



"1 ~"ooJ flop +01 

1602 ri 40^ 



40 

lee^n 40, 

.4 fundamental matrix solution can be chosen with the following behavior at A = 0: 



Vn ^ o. 



*oc/t(A) = [Go + 0(A)] 
Proof: The system (43) is: 

~dA~ ~ 



1 log A 
1 



Go 



1 

9^-9, 2 1 



40oo r x n 



A + A x 



A 



A- 1 



We identify A n + A x and A^ — ^- with £? and Si of proposition 2 in Appendix 1, diagonalizable case 
(98)-(102) (we recall that (103)-(107) never occur when 6^ ^ 0) with a=^ + ^,b = -Of- + fy, 
c = 0. 

The behavior of a fundamental solution is a standard result in the theory of Fuchsian systems. 



The matrix G is defined by G 1 (A + A x ) G = ^ jj q^. 



□ 



Lemma 6 Let r e C. TTie matrices of system (46) are: 
io = 



' , 0o 4 r (r+0 o ) 

' «2 -eg 

_ _ 0o 
' 2 



2 

q2 q2 



/_ r _9o i 4 r (r+0 o ) 



(52) 



There exist a fundamental solution of (46) with the following behavior at /i = oo: 

1' 



*o(/i) 



J + O ( - 



1 log^ 
1 



/i — > oo. 



Proof: To compute ^io an d 4 X for the generic case, we consider the case (105) in Proposition 2, 

0Q I 0x 

2^2 



applied to the system(49). The parameters are a=% + %-,c = #o- In particular, 



Ao-^+A 



0Q + 0x 

2 





1 



(53) 



Here the values of the parameters satisfy the conditions a ^ and a ^ c, namely 8q±9 x ^ 0. From 

the matrices (105), we obtain v4 = £?o + 9o/2 and A x = Bi + 9 x /2. Keeping into account (53), by 
the standard theory of fuchsian systems we have: 



/ 4 () [ - 



+?* ( 1 log/z 
1 



/X — > oo. 



This proves the behavior of 'I'o(m)- □ 

The matching condition ^outW ~ ^o(^) *o (A/a:) becomes: 



ATo(a:) ~ ( e^-ef j_ 

4 fee ri ri 



1 . 0^/1 logx 
1 



From the above result, together with (52), we compute A = K A K 1 , A\ = K AiK 1 . For 
example, 



A = Go 



r + % + ^ logx ^ i og 2 x _ 2 (r + %) logx + 4 



^ log x- (r+%) 



Go" 
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A similar expression holds for A x . The reader can verify that the matching conditions (39), (44) 
are satisfied. 

The leading terms of y(x) are obtained from (50) with matrix entries (Ai)i 2 = — n and: 



r fl2 _ Q 2 



(A))i2 = n 

The result is: 

y(x) ~ a 



5- log 2 x - 2 r + logx 



'o\, , 4r(r + # ) 



2/ 6 fli-flg 



9 2 { 6 \ 

-^log 2 x-2 fr + y J logx + 



4 r(r + 6> ) 



(54) 



4.6 Case 6 ±6 X = 

We consider here the cases (106), (107) of Proposition 2 applied to the system (49). 

Case (106) is the case a = 0, 6o — —0 Xl with a = 0, c = # in the system (49). From Proposition 2 
we immediately have: 

Mf 4)' Mo - 

The behavior of an d ^out, and the matching are the same of subsection 4.5. We obtain the 
same Kq(x). Therefore: 

(A))i2 = n (r - 6» In a;), (Ai) i2 = -n. 

This gives the leading terms: 

y(x) ~ x(r — #o lnx) = x(r + lnx). (55) 

In the same way, we treat the other cases. Case (106) with a = c, is the case <r = 0, 9q = 6 X . As 
above, we find y(x) ~ x(r — # In a;) = x(r — 9 X lnx). Case (107) with a = 0, is the case a = 0, 
#o = — #rr- We find y(x) ~ x(r + O lnx) = x(r — # x lnx). Case (107) with a = c, is the case a = 0, 
#o = 0x- We find y(x) ~ x(r + #o lnx) = x(r + ^ lnx). 

Both (54) and (55) contain more than one term, and in principle only the leading one is certainly 
correct. To prove that they are all correct, we observe that (54) and (55) can be obtained also by 
direct substitution of y(x) = x(Ai + Bi lnx + C\ In 2 x + Di In 3 x + ...) + x 2 (A 2 + B 2 lnx + ...) + ... into 
(PVI). We can recursively determine the coefficients by identifying the same powers of x and lnx. 
As a result we obtain only the five cases (31), which include (54) and (55). 

4.7 No Naive Matching for a = 1 

The condition |Sfter| < 1 suggests that the matching above does not work in the case a = 1 (and 
a = — 1, being equivalent). Let us convince ourselves of this fact by repeating the procedure above. 
A fundamental matrix solution for (40) at A = is non-generic: 

WA) = (G + 0(A)) ( A ; A °.) (J lo f ). (56) 

where: 



G = I I , Vn ^ o. 

40oo ri #oo ri 

A fundamental matrix solution of (46) at \l = oo is non-generic: 

^)) (^O lo SM 1 

where: 

16 r 
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The matching relation: 



K (x) 



(floo+l) 2 -^ 1 / V 

4«oo ri 9oo ri ' 



shows that we cannot eliminate A to obtain Kq{x). 

One case a = 1 is studied in Part III, making use of a non-fuchsian reductions of the system (4). 

4.8 Monodromy Data 

Systems (40) (46) are equivalent to Gauss hyper-geometric equations, as it is explained in Appendix 
1 (make use of the systems (43) and (49) respectively). Therefore, the monodromy can be computed 
in a standard way, using the connection formulae for the hyper-geometric functions. 

We obtain in this way the monodromy of ^out and "J/^. As it is explained in section 2.2, it 
may be necessary to do a transformation ^out ^qut 11 '■— ^OUtCout, in order to match the 
"out" and "in" solutions with a solution * of (4). In this way, the monodromy matrices M , M x , 
Mi of *f> can be obtained. They depend on r. We then compute the traces of MiMj and extract r, 
which is thus obtained as a function of the monodromy data. 

We do not repeat the computations here. One example is the computation of (37) in [17] and [8] 
[9] [11] [5]. 



PART III — Derivations for Non-fuchsian Reduction 

5 Case 9 X = ±0^, 9 X = ±0 O , \mi x ^{A x + A ) = 0. Solution (19) 

We begin by observing that for 6 X = ±# , system (49) may fall in cases (103) and (104). If it is 

so, then A + A x = 0, and therefore A n + A x = 0. More precisely, we start from the following 
hypotheses: 

]im(Ao(x) + A x (x)) = 0, 

x — >0 

Q 

A := lim A x (x) = a constant matrix with eigenvalues ± 

x — 2 

The first hypothesis means that we can write (the trace is zero): 

Aq + A x = [ a . \\ b ^ X } r A , lim a(x) = lim b(x) = lim c(x) = 0, 
\c(x) ± -a(x) J ' x^o v ; x->o y ' x^o v ; 

The second hypotheses implies that the general form of A is: 

/ s + ^ -r \ 
A = ( (s+ejs _ e _e x )> r,seC, r ^ 0. 

\ r 6 2 / 

We also write: 

A x {x) — A=: A x (x), A a + A =: A (x), A + A X = A +A X . 

A x (x) and A (x) are vanishing. We suppose that the slowest vanishing behavior be of order x a °, 
for some a > 0. Namely: 

a(x), b(x), c{x), {A x ) lJ {x), (A )«(a;) = O(x ao ), a Q > 0. 

Finally, we have: 

A 1 (x) = -?fv 3 -(A + A x ) — > -^<7 3 , x^0. 
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5.1 Coalescence of Singularities 
1) THE SYSTEM for ^out- 

We consider system (11), in the domain |A| > \x\ 5ouT . Let us determine the conditions to neglect 
a term x n A x /X n+1 - and all the terms following it - with respect to (A + A x )/X, when x — > 0, 
A ~ x s , 8 < Squt- 



We can neglect 



A" 



< UA q + A x )ij\, Vt,je{l,2}. 



Since \im x ^o(A x )ij are non-zero constants, the above condition is: \x\ n nS <C \x\"° , namely: 5 < 
1 — ao/n. We state this result as a lemma. 



Lemma 7 Let Nqut > 2 be an integer. We can approximate (11) with: 

Nqut-1 



dX 



(A + A x ) Ax 

X X 



+ 



n=l 



A x 
A- 1 



^OUT- 



if and only if: 



Squt < 1 



Co 



No 



UT 



(58) 



Suppose that all term x n A x /X n+ , with n > Nout, have been neglected. We can also make the 
substitution A\ M> —^0-3 in if and only if the error term — A '^Li x 1 is smaller than - "wq !7T j4x ■ 
Namely, if and only if: 



\(A + A x ) i:j \ < 



r-NouT — 1 



(Ac)i 



\Nqut 



This is 1^1^° < | a; |Wo[7T-i-JVo i 7T«oi7T ) nam ely <5 o;7T > 1 - ^2.. 



We can also do the substitution A x n> A, provided that Sour < 1- This is because we can 
neglecting terms X X J+* with respect to A "]| A ' , where both and A + A x arc 0(x CT °), A ~ a; 5 , 
(5 < 1. We summarize the result in the following lemma. 

Lemma 8 Let Nqut >2 be integer. We can approximate (11) with: 



OUT 



dX 



A +A x 



No ut — 1 



+ 



J E (!) - 

n=l 



^00 0"3 

2 A-l 



i/ and only if: 



Nout 

In particular, this means that Squt < 1- 



1 + o-o , 
1 - ~ < dour < 1 



Co 



Nqut 



Example: 

If (To = 1 an d ^V/at = 2 we have: 



dX 



0-3 



x A A Q + A X _6 1 

A 2 A 2 A-l 



(7T, 



< Squt < 2' 



If do = 1 and Nin = 3 we have: 
d^our 



dX 



v x 2 A xA A + A z 
~ + ~ + A 



^00 0-3 



2 A - 1 



^OUT, 



1 x 2 
3 < SOUT < 3. 



2) THE SYSTEM for * /Jv . 
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We consider system (12) in the domain |A| < la;^™. We investigate the condition necessary and 
sufficient to neglect a term X n Ai (and all its next terms) with respect to 4& + xrb- It is convenient 



to write: 

Suppose that A ~ x s , 8 > 8jm- 
We neglect AiX n <= 



+ 



A x A + A x xA Q 



A A — x A — x A(A — x) 





xA Q 


> AiA" 




A(A-x) 


j 


A +A x 
X—x 


> AiA" 



n+2 ' 



\O-0~S vs. I™ 171(5 



Thus, we have the condition 8 > max ^2 1- We have proven the following: 

Lemma 9 Let N IN > 1 be an integer. We approximate (12) with: 



dX 



A A NlN ~ 1 

Ao +^--A 1 ^ 

n=0 



A A — x 



Win- 



if and only if: 



Sin > max 



f <?<> L_ 

IJVjjv + l' TV/at + 2 



We further make the substitution A\ — ^(73. This is possible if and only if two conditions 



are true: 1) 

e 



xAq_ 



X{X — x) 



and 



X—x 
vJV rw -l 



are dominant w.r.t. the term Aq + A x appearing in A\ 



-^CT3 - (A + A x ). 2) lA^^-Mi), i.e. | A^^" -1 ^ | , is dominant w.r.t. the term Aq + A x in 
Esplicitely, the conditions are: 

Ao + A x 



xA 



A(A - x) 



X — x 
> \ A + A, 



^> \A + A x \ (this is always true), 



a;! 1 " 25 > x a \ namely: 8 > 



and 



(JV/ N -1)<5 ^ |™|<to 



> \x\ a " , namely: 8 < 



jA^-Vsl » lAo + A, 
We have: 

Lemma 10 Let iVjjv > 1 be an integer. We approximate (12) with 

d^iN 



~~ 2~ ' 

TO 



N IN - 1 



rfA 



X X-x 2 



A + A x xA 



X-x X(X-x) 2 



71=0 
JVjjv-1 
■(73 ]T A" 



IN 



if and only if: 



max 



1 — (To (To 



2 N IN + 



1' ^VjJV + 2} 



n=0 



< <5/N < 



/AT, 



CO 



As a final simplification, we substitute Aq = —A + Aq 1— )• —A. This is possible if and only if: 



a; A 



A(A-z) 



A + ^ 



A — x 



i +C r -2<5 < la-jTo-^ namely 8 < 1, 



and: 



ir A 



« | A "™-iZ~ CT3 | 



X(X-x) 
We have proven the following: 



1^11+^0-25 c a .(JV JW -i)« i namcly: 5 < 



ao + 1 

^V/JV + 1 ' 
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Lemma 11 Let Njn > 1 be an integer. We can approximate (12) with 

n N IN -l 



IN 



dX 



A + A x xA 6>oo , n 

\(x—x) + — (J " 2. 

v > n=0 



X — x 



IN, 



if and only if: 



max 



1 - (T 



Co 



1 



N IN + V N IN + 2 



< Sin < min 



00 



cq + 1 

Nin-V N IN + 1 



Examples: 

If (Tq = 1 an d Nin = 1, we have: 



JJV 



A + A s 



A — x 



+ 



xA 



A(A - x) 



-C73 



/JV, 



< (5/AT < I- 



If we keep — An instead of A, with no change in the condition on Sin, we can also rewrite: 



in 



dX 



A a 



An 



X X — x 



-CT 3 



IN, 



If (To = 1 and Nin = 2, we have: 
d^iN 



dX 



A + A x xA 6^ ' 

— : 1" T7T r + -^-0-3(1 + A 

A — x A(A - a;) 2 



7AT, 



1 X 2 

3 <<W< 3- 



Equivalently, we can write: 



dA 



A 



A A 



-0-3(1 + A) 



IN- 



5.2 Matching 

We do the matching in the overlapping region \x 

\$OUT < |\| < \ x \SlN_ This 

imposes: Sin < Sout- 

In order for the overlapping region not to be empty, we must choose suitable reductions of (11) and 
(12). If we expect an to be close to 1, we try to match solutions ^out and ^in satisfying one of 
the following sets of systems: 



First choice: 



d^oUT 


' xA 


dX 


/>? 


d^IN 


'An 


dX 


X 



2 A- 1 

0-3(1 + A) 



^OUT, 



The condition to be satisfied for an = 1 is: < 5 IN < Sour < 1 — For (T = 1, this is: 



Second choice: 



d^OUT 


\x 2 A 


xA 


dX 






d^iN 




dX 


A 



A 



2 A - 1 



-(T3 



A A - x 2 

For (To — 1, the condition to be satisfied is: < Sin < Sout < 1 — ■ For <7o = 1, this is: 

1 i / 1 2 

In both cases, the overlapping regions are not empty. The matching procedure will determine 
the leading terms (order x a °) of the unknown matrix elements a(x), b(x), c(x) of An + A x . 
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5.3 Matching for the First Choice: | < 5nv < $out < \- 
We rewrite the systems in a more convenient form: 



A Ju 



dv 
d^iN = 



-xA - 



An 



0, 



CT3 



7 2 z/(z/ - 1) 

2 $oo $oc . Aq A x 

x ir^3 A* + x — CT 3 H 1 - 

2 2 /it jU — 1 

2 $oo ^oo . Ao + A x 

X — CT 3 /•* + X —<T 3 + 



^OUT 
* JJV 

A 



(59) 



/AT- 



, we have v — > oo, 



/x M(M-l) 
Then we substitute A h-> —A in the last term. 

In the matching region \x\~ SlN < \i>\ < \x\~ SoUT , \x\ SouT - 1 < fi < \x 
fj, — > oo. The point at infinity is a non-fuchsian singularity. 5 In order to find the local behavior at 
this point, it is convenient to put the leading term in diagonal form. Let G be the invertible matrix 
such that: 

( 1 

G _1 AG = — -77-03, for example: G 



s+9 x 



and put: 
Then: 

dv 

dfj, 



2 17 00 



^OUT ='■ G ^OUT 

G-^Aq + A^G 6„ 



^out, v -> 00; (60) 



0c 



' ' — CT 3 /" + X — (T 3 + 



A + A 



/i 



(61) 



In order to write the local behavior of ^out and ^in at infinity, we observe that the systems 
(60) and (61) respectively have the following forms: 



dY l 

dz 



' Di D 2 D 3 

z z z z ,} 



Yi, 



dY 2 

dz 



E\ E2 



E 3 



rAz + iAzI — + — + — + ... 



Y2, 



(62) 
(63) 



where and A are diagonal matrices with distinct eigenvalues. In our case: 



^ = X y<7 3 , 

The eigenvalues are distinct iff 6 X ^ 0, Ooc 7^ 0. 

5 System (59) can be also written with ^-03 >-> —Ai: 



A = — a 3 . 



OUT 



dv 



-xA ■ 



A + A a 



-xA - 



A + A x + A- L 



After diagonalization, we get: 



A, 



v(v-l) 

A! 



v - 1 



9, 



dt>oi 



dv 



x-a 3 - 



G- 1 ^ G G- 1 A 1 G 



v-1 



This form is that of a system of isomonodromy deformation for the fifth Painlcvc equation. 
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The theory for such systems is developed in [3] (see also [4]). For any sector of angular width 
7r + e, e > sufficiently small, there exists a unique solution of (62) with asymptotic expansion: 



Yi(z) 



i + — + -4 + ... 



expjfi z} z 



oo. 



fii = diagonal part of D^. 



For any sector of angular width | + e, e > sufficiently small, there exist a unique solution of (63) 
with asymptotic expansion: 



Y 2 (z) 



oxp<j yA z 2 



xA z > z 1 , z — > oo. 



Ai = diagonal part of 



We can always find two solutions Y\(z) and Yi(z) as above, such that the sectors where the 
asymptotic expansions hold are overlapping. We refer the reader to [3] for the general description 
of irregular system with a Stokes phenomenon, and to the Appendix 2 for the computation of the 
matrices Gi, Ki, i = 1, 2, ... 

The systems (60), (61) are isomonodromic. This imposes that Oi and Ai must be independent 
of x. They are: 

Oi = diagonal of (-G _1 (A + A x ) G), 



Ai = diagonal of (A + A x ) 



a{x) 
-a(x) J ' 



We compute: 



G- 1 (A + A X )G=±- 



-(2s + 6 x )a-s{s + e x )b + c 



(-20-56+I) 



7 * V f (2(« + x )a + (s + e x fb - c) (2s + 6 x )a + s(s + 6 x )b - c, 
Since a, b, c vanish, the condition of isomonodromicity implies that: 

Oi=0, Ai=0 

This means that the leading terms of a(x), b(x), c(x) satisfy the conditions: 

a(x) = 0, c(x) = s(s + 9 x )b(x). 

The above conditions mean that if b(x), c(x) — O(x ao ), then a(x) is of higher order, i.e. it vanishes 
faster than x a ° . Note that with this choice of a, b, c we get: 



G- 1 (A + A X )G = 







br 



We are ready to write the behavior of ^out- 
^out = G 



exp ^ x—o?,v f , v — > oo 



G [I + G 1 X + G 2 \ 2 + ...] exp j— <t 3 |J, A^O. 



We use the formulae of Appendix 2 to determine G\. 



(Gi)ij = 2 , J v 

a; ^((ct 3 )m - (o- 3 )jjJ 
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In the second term of the last formula j — 2 if i = 1, j = 1 if i = 2. We compute: 

'-(2a + 6 x ) -2r 



Therefore: 



G-V 3 G = — 

Ox 



(Gi) 



12 



x X 



x "x r 



(Gi)n = -|f (2- + * x ) + fe 2 , (GOa 



-(Gi)n 



(64) 



(65) 



On the other hand, the local behavior of Win is 



O ^ OO oo 

cxp <j a; —0-3 M + x ~2~ cr 3 M 



exp <j — er 3 A + — cr 3 A 



We determine -fTi from the formulas of Appendix 2. 

ifi = diagonal part of (—A) = 



s+% 



— > oo 



A -> 0. 



The matching conditition: 

* 01 /t(A, x) ~ * IN (\, x), x^O, \x\ SouT < |A| < \x\ SlN . 

is restricted to the overlapping sector where both expansions of ^out and ^in hold. Noting that 
^out ~ G and ^im ~ we choose the new solution 'I'oc/t >-> ^out G^ 1 . Then, we expand the 
exponents: 



*oi/t = [/ + GG\G~ 1 \ + GG 2 G~ 1 \ 2 + ...] 



I+-Ga 3 G - + __+... 



The point here is quite delicate. We consider the relation of dominance among terms - and write 
the leading terms of the expansion - as they are in case the G n (x)'s are not divergent when x — > 0. 
Keeping into account that ^fGa^G^ 1 = —A, the dominant terms are: 



*out(X,x) = I + GG 1 G- 1 \ -A% +0 A 2 , 



A 2 ' 

It is important to note that A is dominant w.r.t j, because Sour < \\ namely, A ~ x 5 vanishes 

We expand the exponent in */jv, and keep only the first dominant terms (in the spirit of the 
observation on the dominance relations made above): 

* IN (\,x) =/+%<T 3 X + Ki ^ + o(\ 2 ,^r,x 



A 



^out and "J/jat match in the first term /. We impose the matching of the second term, namely 
the term in A: 

G G x (x) G- 1 ~ ^a 3 , x^O. 
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Namely: 

Gi(x) ~ yG-VsG, x^O (66) 

From the explicit form of G\ and G _1 a 3 G given above, we conclude that the matching is satisfied 
if and only if: 

b(x) ~ —x #00, and a = 1. 

The error in b(x) is of higher order w.r.t. x. The determination of the leading behavior of Aq + A x 
is complete, because c(x) ~ s(s + 9 X ) b(x), x — > 0. Namely: 

c(x) ~ -x {s(s + 9 X ) 6»oo}, a(x) = o(x). 

With such a choice of b(x), one can verify that the terms in ^out and "J/jat which follow the 
second (i.e. which follow the term in A) are actually of higher order in x. Nevertheless, ^out and 
^in match only in the first and second term, being already the off-diagonal entries of the third term 
not matching (i.e. —A and K\ = diagonal part of —A respectively.) 

5.4 Matching for the Second Choice: | < 5 IN < Sour < §• 
We rewrite the systems in the convenient form: 



A X 



<N>QUT 

dv 



-x 2 A v - xA -^±^- e -^^—^ 
v 2 v{y - 1) 



^OUT, 



IN 



@oo , Aq A x 
X ir<T3 + + 7 

2 ii — 1 



IN 



x—a 3 



A + A x 



A 



H /j,(fi - 1) 

Then, substitute A i-> —A in the last term 
We rewrite the systems at infinity 6 : 



IN- 



(67) 



OUT 



dv 



2 6-x O x G {Aq + A X )G 1 

x ttcts v + x — cr 3 — G V 3 G 

2 2 v 2 



IN 



dn 



x—a 3 



Aq + A x 



A | — + — 

^ 2 M 



IN, 



1 1 



Z/, /i — > 00. 



OC/T, 



(68) 



This time the system of ^oc/t is in the form (63), while the system of ^in is in the form (62), 
where: 

Q 

ft = x -^^i ^i = diagonal part of (A + A x ), 

Q 

A = yff 3 , Ai = diagonal part of (-G~ 1 {Aq + A x ) G). 

We impose that ill and Ai do not depend on x, and we get the conditions a — 0, c = s(s + 0^) 6. 
Then, we choose the following solutions: 



^out = G 



oxp J x —a 3 v + x —a 3 



J + ^-G^G- 1 | + GK 1 G~ 1 \ + (j^,x, A 2 ), i/ 



— > oo. 



r Gi 
/+ — + ... 



exp •{ x -^-cr 3 fi \ = I + Gi T + -^-cr 3 A + O ( A ), \i -> oo. 



°The system (67) is in the form of a system of isomonodromy deformation for the fifth Painlevc equation. 
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The relation of dominance among terms are considered as if the G„'s and K n 's do not diverge as 
x — > 0. The matching conditition: 



*oc/t(A, x) ~ *jjv(A, x), x -> 0, 



|<Wt < hi < 



is restricted to the overlapping sector where both expansions of ^out and ^in hold. We note that 
j vanishes slower than A, because Sin > \ (namely, x/X <~ x 1_<5 , A ~ a; 5 , S > 1/2). 'I'/at and ^out 
automatically match in the first term I. We impose the matching of the second leading term, i.e. 
the term in j: 

Gi(x) ~ % G a 3 G- 1 = -A, x^O. 



(69) 



As for Gi, the formulae in Appendix 2 give: 



{Gi)ij - -2 
(Gi)« = -(^)«+2 



(A Q + A x )ij 



i / J- 



X 9oo [{(T3)ii - {<73)jj] 

(Aq + (A a + A x )ji 



X ^cx>[(0"3)jj - {<Ji)ii\ 

In the last formula, j = 2 if i = 1, j = 1 if i = 2. Explicitly: 



(Gi)i 2 = ^ 6, (Gi) 2 i - £i£±M 6, 



(Gi)n = -(A) n - 



b, (Gi) 22 = -(A)22 + 



a(a + Ox) 

•£ Q CyTs 



b 2 . 



(70) 



Therefore (69) 6(x) ~ — x 0oo, and a = 1. As it must be, we get the same result of the 

matching for the first choice. 



5.5 Critical Matching: 



e < 5 IN < Squt < h + e- 



In between the first and the second choice - which hold respectively for | < Sin < <5ot/T < | and 
h < Sin < <5oc/t < I - we can also consider the following approximations of system (4): 



d^our 



dX 



dV 



in 



dX 



xA A + A x 
X 2 X 



Aq | A x 
X X — x 



0-3 



2 A - 1 



^OUT, 



IN- 



Rigorously speaking, the two systems cannot be considered simultaneously when cr = 1. But we 
can consider cr — 1 as a "limit" value - or "critical" value - for the matching of the two above 
systems in the region specified by | — e < Sin < S OUT < | + e, where e > is sufficiently small. We 
write again *o(7T =: G ^>out- Then: 



OC/T 



9x 



G^ 1 (A Q + A X )G 6. 



Y G-V3G 



IJV 



djJL 



&oo , Aq + A x 

x -rr^s H 1- A 

2 



1 

+ 



2+4 + - 



00. 



When we impose isomonodromicity conditions, the diagonal parts of Aq + A x and G _1 (Ao+ Ar)G 
must be independent of x. This gives again a = 0, c = s(s + fl^) 6. Then, we choose the fundamental 
solutions: 



^OUT = G 



1 + 



qOUT 



+ ... 



/ X \ x 
exp 1 x —u-i v \ G = 



= I + GG% UT G- 1 \ + ^GasG- 1 ~r + O (x 2 

2 Ay 



A 2 
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L A 

„<5 



We match them for A ~ x and 8 6 (1/2 — e, 1/2 + e), in the overlapping sector where the 
above expansions hold. Here both A <~ x s and x/X <~ x 1_<5 are the dominant terms. The matching 
conditions are GGf UT G^ 1 ~ % ct 3 and G^ - ^GcrsG" 1 . Namely: 

G? 17 ^) ~ 9 -fG- l a z G, G[ N {x)~-A, x^O (71) 

The matrix k^G^azG can be derived from (64). The matrix G? UT is (65), the matrix G[ N is (70). 
Condition (71) is inclusive of both (66) and (69). Therefore, (71) <=> b(x) ~ — x 9 OQl x — y 0. This 
is again the expected result. 

5.6 Higher Order Terms 

The final result obtained above is: 

-r6oo x\ q 

A + A x = | + o(x), A, = -^a 3 - (A + A x ), (72) 



A * = (s+ejs _ _e )+ oC 1 )' A o = -( (s+ejs _ _ e + °0) ( 73 ) 

\ r S 2 / \ r S 2 / 

Let us substitute the above results into (6). We obtain the first term with no error: 



y{x) ~ . 1 (j , x^O. 

Here, r and s do not appear. Nevertheless, if we substitute in (PVI) the series y = j— g — ^Y^Li b n x n , 
we can compute recursively all the terms, for O = and 9i = ±9^. We find a series: 

1 °° 

y(ar) = - -|ai + Vi„((i;9 M ,«o)/, x -> 0, (74) 

n—0 

where a is an arbitrary parameter. This parameter is actually a function of s, as we prove now. The 
convergence of the Taylor expansion can be proved by a Briot-Bouquet like argument. This will not 
be done here. The reader can find a similar proof in [19] and the general procedure in [14]. 

5.6.1 Determination of a = a(s) 

The system (4) is isomonodromic. This determines the structure A x , A A\ as can be found in [18], 
Appendix C, formulae (C.47), (C.49), (C.51), (C.52), (C.55). If we substitute (74) in the formulae, 
we get a Taylor expansion for the matrix elements, in terms of the parameter a. The leading terms 
have exactly the structure of (72) and (73). We can identify the leading terms to express a as a 
function of s and r. The computations are quite long, so we give the result. When we write the 
leading terms as a function of a and impose that they coincide with (72) and (73), we find: 

_ 9 oa {2s + e x + \) 
a 2(9^ 1) GC 

The higher order terms are Taylor expansions. Explicitly, the first terms are: 

A 1 ^ -^a 3 - (A + A x ) = (75) 



- e -t + (s + 9 x ) S 9 QC x 2 r0 oo {x- V , °° +1) V s+e -- 1) x 2 } 

{(9 x + s)x- (a~-D(»«+;K2«+0«+i)' Y - (** + 



+ 0{x A ) 
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(s + \) - 2(s + e^sO^ x -r{l- (2s + 6 X - 1)6 X x} 
S^±M {i _ ( 2a + 6 X + 1)^ x] - (s + 6 -f) + 2(s + ejsO*, x 



+ 0(x 2 ), (76) 



- (s + %) + 2(s + e^sO^ x r {1 - (2s + 6^6^ x) 

_£(£±M {i _ ( 2s + Bje^ x} (s + %) - 2(a + x )s0oo x 
The above expansions are enough to obtain first two leading terms of (6): 



y{x) = y 



1 eoo(2s + fl x + l) 



+ 



$oo 2(0^ — 1) 

Note that r simplihes. This is the solution (19). 



x +0(x 2 ), 



+ 0(x 2 ). (77) 



5.7 Monodromy Data 

We assume that the matching has been completed as above, and in particular <j = 1. Thus, the 
system (4) can be approximated by: 



d$>0UT 



dX 



A + A x xA x 6>oo 03 



+ 



or 



IN 



A A 2 2 A-l 

A , A x , 0, 



*twr, for |A| > 1x1 s , 5<l: 



dX 



A A 



-03 



* 7Af , for |A| < |af , <J > -; 



The first two leading terms are : 



*OUT = 1+^*3 A + °(f)> A 



- x s -> 0, (5 < 



1 



* /JV = / - A ^ + 0(A), A~x 5 ^0, 5>i-. 

A Z 



(78) 
(79) 

(80) 
(81) 



The above solutions match in the first three terms in the "critical" region A <~ x s , <5 ~ \ (the region 
is restricted to a sector). Namely: 

*oi7t~*jjv~/+^3 *~ A \ + o(x,^,X 2 Y \~x 5 , 5-1 (82) 



Now, for (5 < i, we have: 



A ° + A * ~ x 1 - 5 -». o, X 4^ ~ x 1 - 25 -»■ 0. 



A A 2 
Thus, (78) can be further reduced to: 

oo "3 



In system (79) we rewrite: 



dX 



Av A. 



2 A-l 



A + A x 



^OUT 



xAn 



A A — x A — x A(A — x) 



Then, we observe that, for 5 > \, we have the behaviors: 



A + A x ^ x _ 6 ^ xA 
A — x 



xA 



X(X-x) X(X-x) 



x 1 - 26 -> oo, 



<5> 



(83) 
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Thus, as x — > 0, the system (79) can be further reduced to: 

d^i N xA 



dX A(A - a;) 



IN- 



(84) 



System (83) has the same monodromy of (4) at A = 1, oo. System (84) has the same monodromy 
of (4) at A = 0,x. 

MATCHING # *oc/t = 

We choose ^>out such that it matches with \E" at A — oo. The behavior of this last is: 

1" 



*(A,aO 



I + O 



A 2 



A — >• oo, 



for $ Z. The solution of (83) with the same behavior is: 

*oc/t(A) :=(A-1)-^ 3 



(85) 



As a consequence, the monodromy of \& at A = l,oo coincides with that of iff out- To compute it, 
we consider the loops A — 1 (A — l)e 27 ™ and A i-> Ae . The corresponding monodromy is: 

Mi = Moo = cxp{ 1^6*00 cr 3 }. 



MATCHING * 



OC/T 



OC/T- 



We match ^out with ^oc/t for x/A — > 0. Let us choose the branch (A — 1) = (1 — A) e J 
(1 — A) > for < A < 1. Solution (85) has expansion: 



I+-fv a \ + 0{\ 2 ) 



A -> 



*oc/t(A) = e-''"- 01 
Therefore, for A -> 0, *oc/t matches with *oc/t e -1 * 6 ' 000 ' 3 , where *ot/T is (80). 
MATCHING ^out <-> */jv: This is (82). 

MATCHING * /JV ^ The above matchings imply that e-'f" 00 ' 73 matches with * 

(where is (81)). 

MATCHING */at *jjv: 

In order to determine the monodromy of system (4) at A = 0, x, we need to find a fundamental 
solution of (84), that matches with ^ JJV e~ % ^" 3 for x/X — > 0, where VP/at is (81). A fundamental 
solution satisfying these requirements is: 



*/at(A,x) := G 

Actually, this has the behavior: 
# /Ar (A,x) 



('-!)' 



Q-~L g-if SooO-3 



where ^f-Ga^G 1 = —A. The first two terms match with e l 2 



i£ #00 CT3 



A^°' 



, as required. 



As a consequence of the matching, the monodromy matrices of "J and at A = 0, x coincide. 
To compute them, we write the local behaviors (for i^O fixed) of Win'- 



$ IN (\,x) = 



Gx^ 3 



[l + 0(A-x)] (A-i)-Vs [ G -i e -ifeoo<r 3 ] ; A — > 



G(-x)-^ L<J3 [1 + 0(A)] A^ ff3 [G-i e -*f0oo<x3] ) A ^ . 
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It is not necessary to specify the branch of (±x)~~5" J3 . We resctict to the case 9q and 9 X Z, so 
that the matrices i?o, Rx in (17) are zero and the matching of ^ with the above behaviors of ^>in 
is realized. 

We can compute the monodromy matrices for A i— > Ae 2 '™ and (A — x) h-> (A — x)e 2 '"' 1 respectively: 
Mo = e^ e ~ CT3 G exp{i7r^<T 3 } G^e-^^ 3 , 

M x = e ^ 9 ~ CT3 G exp{-i7r6> x a 3 } G^e-'i ^ 3 , 

Note that Mi = M^ is invariant for Mi i-» e _i ^ ^^Mie'^ e ~ <T3 . With this fact in mind, we obtain 
the result of theorem 3, point b). In particular, computing the trace of Mo Mi, we get: 

9 X [2 cosM^ + X ) ) - tr(MiMo)] 

s = 



2 [003(^(6*00 - e x )) - cos^^ + e x ))] ' 
6 Case a = ±{8 1 -#oo), ±(0i + 0oo). Solution (18). 

This case shows new features, namely r and r\ may be functions of a;. For a = ±(6>i— #00), ±(6 , i+6 , 0O ) 
the matrices (41) (42) become: 

a = ±(9 1 -9 oc ): ii=(~^ I 1 ), A + A x =(^f* (86) 
a = ±(^+000): A 1 = M A + ^=( Q 2 ^ 



2 / \ " 2 

r 2 /q \2-ir 2 /q ,q \2j 

The transpose matrices may be considered (namely, re-define n >->• — -1 ^ — n )■ 

The matrices A and are again (47) and (48). For definiteness, we will consider the upper 
triangular matrices A\ and A a + A x , and the choice a = 9\ — 9^. We distinguish three cases: 

(I) n is a non zero constant: In this case we just repeat the general matching procedure and 
find y(x) as in (51). r 7^ is constant. 

(II) n = and r constant. This is a very easy case, because A\ = — ^03, A + A x = 03. 
Therefore a fundamental solution of (40) is: 

01 fit ~ 1 01 -0QC 

A ^ CT3 , A^O. 

The solution of (46), with r constant, is the same of the general case. The matching is possible as 

" -flop 

2 

x 2 



*oc/t(A) = (A - l)-^3 A ^-^3 (_!)-^3 



7+ya 3 + 0(A 2 ) 



in the general case, with = I X 2 OoP-oi I . Therefore 



16<7 2 r X 4a 



We obtain: y(x) ~ r/ - - — ^"l 12 , , 5 ; = , x L-7 ",Z n = 1- This is the singular solution 

yy > s[(Ao)i2+(Ai)i 2 ]-(Ai)i 2 x[ x»i «°o +0 ]-0 y 

y(x) = 1. 

The case r\ — r x (x) — > may give a non-singular j/(x), provided that also r = r(x) — > 0. This 
will be proved below. 

(Ill) Case ri = ri(a;). 

A priori, we don't know if it is possible to have an x-depcndent r\. This is not actually possible 
if we approximate (4) with the fuchsian system (40) (i.e. (13)). Actually, (40) is a system with 
reducible monodromy (upper triangular): 

-^r = \{ h + % a-t w 
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There is a fundamental solution (obtained by variation of parameters): 



^out = 



(1-A)- 




(1-A)- 



^ riF(l-0i, 000-01,000-01 + 1^) 

&oo— 01 

1 



«1 -»o 







A 





. »1- 



Here, denotes a Gauss hyper-geometric function. This solution has a diagonal monondromy 

matrix at A = and an upper triangular monodromy matrix at A = 1, with n in the (1, 2) matrix 
element. Therefore, r\ must be independent of the monodromy preserving deformation x. We are 
sent back to the case (I) of constant n . 

The only possibility for r\ to depend on x is that the matrices Aq, A x have a behavior, for x — > 0, 
such that the system (4) is approximated by a system (11) with singularity of the second kind at 
A = 0. Namely, system (4) must be approximated at least by: 



d&our 

dx 



xA x A +A x A x 
A 2 A A - 1 



OUT 



(87) 



Hypothesis: We consider the case: 



xA x — ^ + higher orders x — > 0. 



In the above hypothesis, (87) is: 

d^ouT _ ( ( p(a;) \ _j_ 
dA ~ 1 1 / A 2 



+ 



2 




l-0o 



+ 



A- l 



(88) 



Let us write a ^out as a vector ^ ^ J . The system becomes: 

dtp! (0!-6 oo 0! \ /n p n \ . dy> 2 /< 



'i 



2A 2(A-1) 



y>2. 



This system is solvable by variation of parameters. Let C\ and C 2 be integration constants. The 
general solution is: 

^a(A) = C 2 X^ X {X- 1)^, (89) 

P F(i-e 1 ,e oo -e 1 -i,0 oo -e 1 :X) 



ipi(X) = a A^^(A - 1)"T + C 2 e MWl 
n -p 



+ 



o-01-l 

F(l — 0\, 0oo — 0\, 0oo — 0\ + 1; A) 



A 



A - 3— (A - 1)- 



+ 



(90) 



Here F(...) denotes a Gauss Hyper-geometric equation. The choice of the branch is such that for 
< A < 1, we have < 1 - A = e" <7r (A - 1). 

In order to write the local behavior for A — > 0, we expand the hypergemoetric functions and 
(A-l) ± ^/ 2 = (e-(l-A)) ±ei/2 : 

i I + 

1 + 0\ — #oG A 



^ = Ci e^A^^ ( 1 + yA + 0(A 2 ) ) + C 2 e l ^X^ 



P 0i(0i-0 OD + 2)-2r 1 (0X-0QQ + 1) 

2(0i-0 oo )(0 1 -0 oo + l) + 1 J . 



^ 2 = C 2 e*' 81 ( l-yA + 0(A 2 ) ) A" 
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Therefore, we can take as a fundamental solution the following matrix (choose C\ = e 1 ?® 1 , C 2 
e -ifOi) : 



^OUT = 



(1 + ^A + ...)A 



p 1 , P 8i(9i-9 30 +2)-2n (81-8^ + 1) 

l + fll-floo A "I" 2(01-0oo)(01-0oo+l) 



#00 — #1 

A 2 



(l-fA + ...)A- 



n P \ /-, p 9i(9i-9oo+2)-2ri (9i-9oo + l) 

l+ej-floo \ l / 1 2 (e 1 -e 00 )(e 1 -e oc +i) 

o o / A In i 



O(A) ^ A— 



6.1 Matching 

The above solution must be matched with the solution of system (46), with a = 6\ — 9 Q 



e 1 -e a 

2V '(£)^~ 



00-81 
A\ 2 



From a standard computation we find: 

/ flp — (S x — a) 2 



4<T 



<r+l 



16cr 2 (cr-l)r 4<t(<t+1) 



(9 Mfl*-^r)(eo-(9*+^r) (8x+^) 2 -9 2 



Note that r in (48) is constant, because the monodromy of (46) depends on r. 
The matching relation ^out(^) ~ Kq(x)^o (A/ar), reads: 



n e 

i+ei-6oo \ l 



i P 9i(9i-9 3 o+2)-2r 1 (81-800 + !) 
2(6>l-0oo)(01-0oo+l) 



+ 0(A) 



# (a;) 



/ + #i- + ... 



x 2 

a: - ^ 





81 -»o 



Namely: 



P 8 1 (8 1 -8 oc +2)-2r 1 (81-8^+1) 
2(8i-0 oo )(8i-0 oo +l) 

1 



i^o(^) 



x 2 

9l - 9o 
X 2 







l+8i-8o 




The first equation above is: 



xK (x)F^ A 
(J x 2 



ifo(ar) 



1 







P 8i(8 1 -8 ao +2)-2r 1 (81-830 + !) 
2(fli-e oo )(8i-e oo + l) 

1 



a; 2 



We substitute this in the second equation. For simplicity, denote the matrix elements of K\. 
We obtain: 







I+81- 





where: 



(*) := 



x K n + (*) JT 21 a;^" ^ 1 K 12 x ^-^ +1 + (*) a; if 22 
#21 a; -- ^ 1 x if 22 

pgi(gi-goo+2)-2ri (gi-gpo + 1) 
2(^1-^00)^1-^ + 1) 
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The element (2,1) must vanish. This occurs iff: 

x 0oo-0i+i 0j forx^O; <=^ - 6»i) > -1. (91) 



We substitute this result in the element (1,1) and then we impose that it vanishes: 

i - floo + 2) - 2ri (gi - C 

2(01-000X01-000 + 1) 



P 0i(0i-0oc + 2) -2n (01-0QO + 1) K ^ x e^-e 1+ i ^ Q 



This implies that p 9 1 (9 1 -6^+2)- 1r x (9 1 -6^ + 1) = oix ^ 6 ™- 1 ). 
From the clement (1, 2) we have: 



#12 x e ^ +1 + o(x e ^) K 22 . (92) 



l + 01 — 0QO 

This relation may be satisfied in two ways: the first is that p = p(x) = o(x 6l ~ eoa ). The second is: 

p = p(x)~-(l + 8 1 -e oo )K 12 x e i- e -+ 1 =-rx e i- e °° +1 . (93) 

In both cases p is a function of x. We are going to prove that the monodromy of ^out is independent 
of p(x) (namely, of x) if and only if: 

01 - 0OO + 1 , n ,N 

p = n . (94) 



This fact rules out the first possibility, because (92) becomes: 

#12 x 6l ~ 0ac+1 + constant x p(x) x, 



P{ x ) ,-, II, II . + ! 



0oo — 01 — 1 

so the last term in the r.h.s. is a higher order correction and p is given by (93). Before proving (94), 
we complete the matching procedure. Using (94), we find: 

A " ( ^lo . j 1 o ^y 

We are ready to compute A = K AqKq~ 1 , where A is (47), for a = 6\ — 0^: 

/ (io)n + g(A ) 21 x - 2g(A ) n - g 2 (A ) 21 x) x e ^+ l \ 

= V J . (95) 

V (Uo) 21 -(io)u - 5(io) 2 i x J 

The first term of each matrix element certainly contains no error. 

We can now substitute (i )i2 = (-£ + 0(1)) x 6l - e °° +1 and A 1 = -r±{x) ~ fll _% 1 3o+ i r:E 9l ~ e °° +1 
into (6), and find: 

/ N 01 - 0OO + 1 „ 

y( x ) ; a ' x ~* °' 

t — Poo 

Proof of (94)- We compute the monodromy of the solution ^out- At A = this is given by the 
matrix Mq UT = exp{i7r(0i — 0oo)o"3}, obtained by the expansion of ^out at A = as we did above, 
with the choice d = e^ 01 = 1/C 2 . 

Let us study the monodromy at A = 1. We need to expand (89) and (90) at A = 1. First of all, 
we use the contiguity relation: 

F(l - fli, Ac - fli, floo - fli + 1; A) = ^ I [(1-A)" 1 -F(l-fli,floo-fli-l,floo-fli;A)] . 

This is used to rewrite tpi: 

P-n (l-X) 6 



i>i = CiA^^(A - I)"- + C 2 e 17101 A^^(A - 1) 



0oo — 1 A 
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+ 



n 15 7 + t) T F(l - Oi, Ooo - 0\ - 1, - Oi; A) 

Coo — i>\ — L Wqo — 1 / A 



Then, we substitute in f/'i the following connection formula: 

F(1-6 1 ,9 QC - 6 1 - 1,6 00 - 9 1 ;\)= r( ^ (g °° ~ 6l) F(l - 1; 0^ - 6 1 - 1, 1 - i; 1 - A)+ 

Thus, tpi has the following structure, when A — V 1: 



V>i =Cie- 4 2 & i A - 5— (I-A)-t + C<2.e l * A" 



(1-A)^£a„(l-A)"- 



n=0 



where a„,6„ ^ are the coefficients that follow from the expansion 1/A and the hyper-geometric 

_l_ Coo -Hi 

functions at A = 1. When A* 2 is also expanded at A = 1, we find that ^out has the following 
structure: 



(seriesx) (1 - A)"* (series 2 ) (1 - A)* + (j^-, + (series 3 ) (1 - A)"* 



*&OUT = 



(series4) (1 — A) 2 



Here "series" means a series of the form X^^Lo C ™(1 — ^)™' w ^ n c o 7^ 0. We just give the dominant 
term, with the choice C\ — e l ^ dl — I/C2: 



*ot/r = [I + 0{1-X)]x 



rX-MT(g°o-gi) 



where: 



( X A ) 2 _fc^T + (**)r(i-ei)r(e 00 -ei-i) 


(**) 



(i-A)^N«?(i 

(1-A)* 
P ri-p 



do — #1 — 1 $co — 1 

Let be non zero arbitrary numbers. We conclude that: 



^OUT = 



X 



w 
(l-A)-^ 



T + (**)■ 



r(-ei)r(e 



[7 + 0(1 -A)] x 







'r(i-0i)r(0oo-0i-i) 

V V K**) ] 

n u 2 ^ I \ Q 1 \ p- r ^ 1 r-i-i-i r (- 9 i) r (^- fl i) 

U (1-Aj2yyu w ^ oo _ 1 + (**) r{1 _g l)r{ g oo _ ei _ 1) 

The monodromy matrix is then: 



r^-i) 



Mf* 77 = Q-^-'^^Q, Q := 



1 i^IMfcM 

« ^* ' r(9oo-i) 







p—ri 



e—^i + (**)r(i-ei)r(6oo-ei-i) 
This is independent of if and only if (**) = 0. This proves (94). □ 
(94) simplifies considerably the structure of (90): 

^ = d A^ 2 — (A - I)-- - C 2 1 (A - 1)-. 

"1 



r(-e 1 )r(e 00 -e 1 ) 



3G 



With the choice C\ — exp{z-|#i} = I/C2, we finally h 



we: 



^OUT — 



A^(l-A)^ 



JVoie: The monodromy at A = 0, 1 is diagonal: M§ UT = cxp {in(6i - 6*00)0-3}, Mf 3 ^ = exp {-iwO^}. 
It is independent of r\. In case (I) - namely, r\ is a non-zero constant and the system for is 
fuchsian - if Mq UT is in diagonal form, then M° UT is upper triangular and depends on r\. 

6.2 Higher Orders Terms 

We may repeat the same procedure of section 5.6. We write a Taylor expansion y(x) = (6\ — 9^ + 
1)/(1 — 6*00) + X)n>i bnX n , substitute it into (PVI) and determine recursively the &„'s. Then, we 
may substitute the result in the matrix elements of Ai, i = 0, x, 1, according to the formulae of [18], 
and find the higher orders of the matrix elements. We just give the result. 

1 a \ 1 A \ (floo - 0i) 2 + Oj - Op 
(^)n = - (A x ) 22 = — - ^ + 

, 6>i [(gi - goo) 2 - (go - e x ?Wi - M 2 - (8 + e x f] 2 
+ 8 (ft - M 2 [(0i - M 2 - 1] I+ 

s fl giKgi - goo + 2)((9i -gpo) +gg -0 2 ] ) ei-Soo+i 

( ^ )l2 " 2(0 1 -0 OO )(0 OO -0 1 -1) + * 

. 1 [(0i - e x ) 2 - (6 Q - e^w-, - e^f - (go + 0^] r 

( ^ 21 = r 16(0oo-0i) 2 I + 

_ 9i[{9qo — 0i){0oo — gi + 2) + g| — 2 ] ^ 2 _j_ Q^ a .3> 



2(6*00 - 0i ) 3 (e oo - 0i + 1) v '] x«i- e -+i 

^ 4 \ _ (a \ _ (6*00 - 6>i) 2 + gj - 2 
(A )n - - (Ao) 22 - m _ 6ac) + 

0i [(81 - 0oo) 2 - (go - :c ) 2 ][(0i - M 2 - (go + ^) 2 ] , n( 2, 
8 (0i-0oo) 2 [(0i-0oo) 2 -l] * + 

• < 1 . ( 1 6>i[(6>i-6> 00 ) 2 + 6> 2 -6> 2 ] n( ) \ fli-Soo+i 
(A))l2 - r i x 2(0i - M(0oo - 6 1 ! - 1) + ( } f 



i [(0i-0oo) 2 -(0o-0.) 2 ][(0i- M 2 - (ft) + ^) 2 ] 

r 16(6»oo-6»i) 

\2 1 /)2 /)2i 



(^0)21 = — Q ^ 1 X + 

01 [(01 — #oo) 2 + 8'o — 2 ] ■> , . ; 



2(0oc -0l) 3 (0oc -01 + 1) X +0(X ' ) j a;fll-9oo + l 

V _ M x _ 8l 6'l[(01-0oo) 2 -(0 t ,-0 :c ) 2 ][(01-0oo) 2 -(0o+0 :c ) 2 ] a 3 
(Al)n - -(A l)2 2 —- 16[(0oo-0l) 2 - l](0oc -0!) 2 X + ] ' 



1 lj " " 0oo-01-l\ + 2(0oo-0l)(0oo-01-2) I+U(1, J" 



-foo+1 



(Ai) a i 



0l[(01-0oo) 2 -(0()-0a ; ) 2 ][(01-0oo) 2 -(0O+6' a; ) 2 ] 2 \ 1 



16r(0oo -0i + l)(0oo -0i) 2 J x^i-^+i 



The above leading terms of y(x) are related to the above formulae through (6). The truncation of 
(Ao)i2 and (Ai)i2 above is enough to reproduce the first two terms of solution (18): 

/ , 6*i — + 1 Oi[(0i - 6>oq) 2 + 6l — 0q + 26 1 - 26>oq] 2 
2(1-^)^-^)^1-^ + 2) 
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6.3 Transpose Case and General Result 

In the above computations, the condition x e °°~ ei+1 — > was necessary to do the matching. We can 
repeat the matching procedure starting from the transpose matrices: 

A '=(~t £)■ -40' x *•-{*■) o)- 

The procedure is exactly the same, with the necessary condition: 



x 



01-0oo + l 



— > 0, when x — > 0. 



As a result we obtain again exactly the matrices A , A x , A\ above. The reader can convince himself 
without doing any computation, simply looking at the structure of the first terms of above matrices. 
For example, let us have a look at A x . Denote the constant terms with letters ci, C2, etc. We 

ci + ... + ^i-flcc+i' 

|c3 X + ■ ■■} ^-L + i -Cl + ... 

Cl + ... { c 2 » + ..■} ^-^ + 1 

{a + ...} a;9«,-9l+l -ci + ... 

The role of 6^ — 6\ + 1 and #i — + 1 is just exchanged. By continuity, the matrices A n , A x , A\ 
computed above hold for any value of 9\ — 9^ £ Z, 9^ ^ 1. 

The matching procedure can be repeated in the same way in case a = — (9^ — 9\) and in the 
case o = ±(#00 + #1), which yields (21). For this last cases the results are obtained just by the 
substitution 9\ i-> —9\. 

6.4 Monodromy Data 

We compute the monodromy data for the case a — 9\ — $007 the other cases being analogous. In 
this case, the matching has been realized by: 



* ol / T (A,x) := 

A"^(l-A) ! 




*itr(\x) = K (x)^o [-), * (/i) 



/ - () [ - 



«1 -Baa 

-<T 3 



/J, 2 % jU -)■ OO. 



Let *(A) denote the solution of the system (4) of (PVI), such that: *(A) = [I + 0(A -1 )] A--? L,T3 , 
A -> 00, ^ Z. 

MATCHING * «• ^out- 

With the choice 1 - A = e~ iw (X - 1) (1 - A > for < A < 1) we have: 



C/T 



I -Oil 



\-^^ e J t Wl<T3 , A ^00 



A / 

Therefore, the correct choice for $out, which matches with ^ is: 

^our h '■— ^outCout, Cout ■= e ^ 0icr3 . 
As a consequence we obtain: 

i) the monodromy of "J/ at A = 1, 00, which is equal to the monodromy of ^Q^ h : 

Mi = e'^ 01 " 3 , Moo = e^ 9 ^ 3 ; 
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ii) the correct choice for ^f^ tch ■= K (x) * Cqut- 



MATCHING # <-> qiMatch 

This is realized by construction. A consequence, we can compute the monodromy of ^ at 
A = 0, x. In order to do this, we need the local behavior of ^fj^ tch (X, x) at A = 0, x. We start with 
recalling that: 

* (ao = ^ - i)**ooo, <w - 

Here, <pi, ^2 are two independent solutions of a Gauss hyper-geometric equation (see (108) in 
Appendix 1): 



d<p 



fi(l - fi) — \ + (1 + c - (a + [b + 1] + 1) n) -f - a(b + 1) tp = 0, 
(III- '" 



where a := f + ^ + 2f - 6+ 1 := ^ + ^ + ^ - 2f + 1, c+ 1 := O + 1- The functions £1 and 
£2 are obtained from <pi and 932 by: 



r 



dp 



(96) 



In order to have generic solutions (i.e. non-logarithmic solutions) of the hyper-geometric equation, 
we must require: 

Then, we have the following sets of independent solutions at fj, = 0, 1, 00 respectively (we denote by 
F the Gauss hyper-geometric function): 



(o) 
<Pi 



(0) 
¥>2 



= F(^ + ?f + if-^, %^ + %- + %^ — + 1 + O ; m) , 

M r \ 2 2 ' 2 2 ' 2 2^2 2' 1 ' 1 u 0, ■ 



4 1 )=F(% + % + 



£1 Z£L 

2 ' 2 



#1 #t> 



+ 1, 1 + X ; 



¥>2 



(i) 



(00) 



2 2 2^2 



(00) 

<^2 



fpffoj-ikj-foo-iL <k _ <>o 4- £0° _ #1 a _a.i\ 
" J ^ 2 ~ l ~ 2 " r 2 2' 2 2^2 2 ' 1700 171 ' ^/ ' 

The connection formulae can be found in any book on special functions: 

Cbi, - 7r < arg (1 - n) < n. 



(0) (0) 



<Pl ' ^2 



(oo) (oo) 
^2 



C oo, < arg fi < 2n. 
Here the connection matrices Coi, Cooo are (34) and (33) respectively. 

From the Taylor expansion of the hyper-geometric functions in <p|°°' ) and (96) we compute: 



(oo) 



l + O 

I + I - 



// 



— a - 



/i 



l + O 



/i 



A* 2 (M-l) 2 C 2 = 



l + O - 



// 
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It follows that the matrix ^o(fi) = 1)~<3>o(m) with the prescribed behavior 7 + [-) 



#1 — #oo _ 
/-i 2 



at /i = oo is: 

(oo) (oo) 



4i 4 2 / V u fli-0oo+i 



*0(M)= M 2 (M-l) 2 (J(co) g») J ; 

Let: 

V u ei-eoo+i/ V u ei-Soo+i/ 

We conclude that: 

*jr*(A,*)=^o(*) (-) f--0 u " i • 

\* /^ 0) (i) ^ 0) (i)', , 

= K (x) - --1 CSiC, 

= K Q (x) - - - 1 CoiC -^C. 

The behaviors of the above matrix at A = x, is easily computed from: 

fa, (i) 
H 2 (m- 1) 2 v>i - 



A\~ /A 



yu 2 (fj, - 1) 2 1 - fi) = x 2 2 FI...; I A 2 (A - x) 2 , 



fa, 1X £» (l) 
/Li 2 (// - 1) 2 <p\> 



fa, / -?0 + ix / x-A\ .fo 8* 

At 2 (/z - 1) 2 (1 - /i) 2 = (-1) x x 2 + 2 F ...; A 2 (A — a;) 2 ; 



and: 



fa, 1N »2L (o) 
M 2 (/" - 1) 2 ¥>i 



fa, .x^r-,, x _^i_hL „ / AA >£*.fa 
= M 2 (A* — 1) 2 ^(-;aO =x 2 2 F (A-x) 2 A 2 . 



fa, ..Sjl (o) 
At 2 (m - 1) 2 = 



= - = x^-^ F ( ...; ^ J (A - x) + \-*. 
From these we compute: 

K (x)^(p-1)^ M°j ^)=4 N (x)(I + 0(\))\^, A^O 

V SI S2 / 

[i~ [p. - 1)* % % = V^W(/ + 0(A - x)) (A - x)^, A -> x. 

A SI S2 / 

Here, we don't need to explicitly give the invertible matrices iI>q N (x) and ipl N (x). From the above 
procedure, we find: 



M = CT 1 (Cooo exp{i7T0o<73} C^) C, M x - C" 1 [c 0o o ( Cof exp^^as} On) 



C. 



We finally observe that Mi and are diagonal, so they are invariant for the conjugation M i-> 
CMC" 1 . With this in mind, we get the result of theorem 3, point a). 

We stress that A + A x and A\ are upper (or lower) triangular matrices, and the group generated 
by M X M and Mi is reducible. 
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7 Case 9^ = 1, 9 X = 0, a = ±1. Solution (20). 



A + = 



In section 6 we imposed that 8^ ^ 1 and O^—Ox £ Z. Here we consider the case ^ = 1, 9^— 9\ = 1. 
We have: 

' ~\ r A a -f° ~ ri 

\ ) ' Al - ^0 

Also the transpose matrices are possible. For |A| < | a? | 5jf " , we use the reduction (46) (namely, (14)). 

The matrices A , A x are given by (47) and (48) with the substitution a = 1 or —1. For definiteness, 
let us take a = 1 in the following. 



For |A| > \x\ SouT , we approximate (4) with: d ^°» T = 



xA x , A +A x 



A-l 



^oc/t- For defi- 



niteness, let us consider the case when A\ and A + A x are upper triangular. Again, we make the 
hypothesis that the leading terms in xA x define an upper triangular matrix: 



p 




higher orders. 



Therefore, we will study: 



OUT 



A 2 V / A 



-2 7 

i 



+ 



1 







-n 



A- 1 V 



OC/T 



This is a reducible system. To solve it, we write ^out m vector notation ^ot/T 
system becomes: 



The 



Ha~ = ~2A 1 A 2 A~ ~ A - ~T 



dih 1 , 
dA = 2A^ 2 - 



^ 2 = C 2 A^, Ci,C 2 eC. 



The solution obtained by variation of parameters is: 

^ = CiA-5 + C* 2 A-5(pln A - n ln(A - 1)), 
We can choose the following fundamental matrix: 

/\-\ A-5(plnA-nln(A- 1)) \ 

V A^ J 

Its monodromy relative to the loops A n> Ae 27 ™ and (A — 1) (A — l)e 27 ™ is respectively: 



\~2 







1 

A~ 



1 p InA 
1 



ri ln(A - 1) 
1 



-1 





-27rip\ 



^OUT ^OUT 



-2niri \ 







J- 



Therefore, p and n mwsi be independent of (the monodromy preserving deformation) x. We observe 
that any fundamental matrix solution of the form: 



^ OUT 



1 /(*) 







1 



has the same monodromy of ^out, for any arbitrary function of x. This fact will be used soon, 
with the choice f(x) = —p\nx. 

With the choice of the branch of ln(A - 1) = ln(l - A) + iir (i.e. (A - 1) = e 47r (l - A), 1 - A > 
for < A < 1), it is convenient to redefine ^out by: 



^OUT '■- 



\~2 

A^ 



1 p InA 
1 



1 -^ln(l-A) 
1 



1 


A~5 





-n ln(A - 1) 
1 



1 iirri 
1 



Therefore: 

^OUT 



1 n 

1 



+ 



i7rri 




E 

71=1 



A n 
n+ 1 





A* 



A-5 




1 p InA 
1 





A^ 



1 p In A 
1 



A 
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7.1 Matching 

The solution (^) has been introduced in section 4.7. 



I + O 



o 1 



Some adjustments are necessary. Let us consider the permutation matrix P := ^ ^ j and redefine: 

*o^*oP(j /i), 9out»9out(1 ~ P ^ X ) (j p °i). 

As we have already observed, this re-definition does not affect the monodromy of ^out, which is 
independent of x. The matching relation ^outW ~ K (x)^o (V x ) becomes: 

1 n\ fl 0\ ( \-\ \ (l In 



1A0 i.M a*Mo 



The matching is thus realized, with the choice: 



1 OWxi \ / A-2 \ / 1 In 4 



^ M rUl A3 M 1 



It follows that: 

-■ 2 -flg-l ri [l-(9 -8 x ) 2 ][l-(8 +9 x ) 2 ] p , „ >j [l-(g -e x ) 2 )][l-(fl +9 x ) 2 ] _ ' 

16 x x^' 1 2 p 16 X 



A 



A x 



1 [l-(9o-ex) 2 )][l-(9o+M 2 ] r 2 -fl 2 -l ri [l-(9 -e ;c ) 2 ][l-(9o+9 ;c ) 2 ] 

p 16 ^ 4 p 16 ^ 

e 2 -e 2 -i D [i-(e -e ;c ) 2 ][i-(go+8 :c ) 2 ] T . p , „ el-e 2 Q +i r \ [i-(e -e x) 2 ][i-(e 0+ e x ) 2 ] 

4 p 16 X a; "•" ' 1 2 ' p 16 ,L 

1 [l-(^o-gx) 2 ][l-(9o+gx) 2 ] „ gg^gg — 1 ri [l-(g -g x ) 2 ][l-(g 0+ g x ) 2 ] 

p 16 x 4 " r p 16 ^ 



The leading term of each matrix element certainly contains no error. Note that r has simplified. 
Actually, the constant p plays the role of r. 

We will not repeat again the discussion for the higher order terms. The matrix elements are 
Taylor expansions, corresponding to a Taylor expanded y(x), the convergence of which is proved by 
a Briot-Bouquet argument. The first two leading terms of A and A x above are actually correct. In 
particular, we need: 

(A))i2 = -f +r /° 2 ~^ + 1 +0(x). 

The leading terms of A\ are: 



r^ li-(e -e x m-(e 0+ e x f ] x 2 + o{x3) _ n + r ^l^l x + 0(a . 2) 

ri [l-(9o^e x ) 2 ] 2 [l 4 -(e +9 x ) 2 ] 2 3 ,4 , Q( yS) n [l-(eo-e x ) 2 ][l-(eo+9 x ) 2 ] ^2 , Q ^3) 



The above truncations of (Ai)i 2 and (A )i2 corresponding to the first two terms of the Taylor 
expansion of the solution (20), through (6): 

y(x) = a + ^—-^{l + el-6l)x + 0{x 2 ), where a := ( 1 - — \ . (97) 



P 



Observe that a depends on the monodromy datum 
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7.2 Monodromy 

The matching ^P/jv <-> ^out has been realized by: 

\~i 0\fl p -n ln(A - 1) \ / 1 -Inx+in^- 



^OUT = 



and: 



A3 J \ 1 

R- 1 



1 







R 



1 







[(e + e x ) 2 - i][(e - e x f - 1} 

16r 



*o(M) = 



I + O 



Namely: 



^in ~ K (x)P 



X 2 



iT 1 a;- 2 



for = A/x ->• oo, |A| < |x| 5 ' N , x ->• 0. 



/i2 





A"5 
As 



1 a 

p 



1 

j? in/x 1 r 



1 ln^ 
1 



A"3 



/i — > oo. 



1 lnf 
A* M 1 



MATCHING * «• *oc/t- 

The correct choice of "J/qc/t must match with: 



A-T^ A i ; A ^ oo; 0, 



oo 1 7 



where: L = ( jj ^ ri ) . This form of L follows from the standard theory of Fuchsian systems, 



and from the expansion of the system (4) at v := j- — > 0: 



dv 



~ - {A l+ xA x ) + 0{v) 



v ->• oo; 



Thus L12 = — (A\ + xA x )i t2 \x=o = r\ — p (this is computed from the expansions of Ai and A x at 
x = obtained before). We expand "fot/T at A = oo. We easily get: 

i i (\ (o-rA lnA\ /I JTT^-lnx 

*OUT = lI + O(\- 1 )]\-^^ (P 7 lllA j % 

Therefore, the correct choice is: 



*OUT ■- *0 



UT 







A _ fCT3 , 1 plnA^ fl -nln(A-l) 











As a consequence, we obtain: 

i) the monodromy of * at A = 1, oo. This coincides with that of which is easily computed 

from the local behavior: 



*#$? fc (A) = [J + 0(A-l)] 



1 -nln(A-l) 
1 



A4 1. 



Thus, for (A — 1) i— >• (A — 1) cxp{27ri} and A i->- Acxp{27ri} the monodromy is: 

1 — 2niri \ Jjf ^ — 1 2i^i{r\ — p) 



Mi = 







1 



M 







-1 
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ii) the correct choice of ^in, which matches with "J. This is: 



1 17T— — In x \ 1 f\ n l n r — iirn 

p , = * /A r(A,x) f P lnX lnn 



1 I -I« V ",«V I Q p 



K (x)^ (p) Cin, C IN 



p/R 

1 pin a; — inri J ' 



Important Remark: Cjn depends on x. But the monodromy must be independent of x. For this 
reason, we will have to consider the substitution 

yMatch yMatch( CiN yl = Ka{x) y o{p) . 

This makes the monodromy at A = 0,x independent of x. The corresponding transformation: 

,t, Match , . *j,Matchf-i — 1 
^OUT ^ ^OUT C IN , 

changes Mi, M x , but does not introduce a dependence on x. Namely: 

IVh ^ CiNMiCiN^ 1 = { „ L „ ?), Moo m- C/jvMooC/jv -1 ' 1 



2^11 i? i^. iH ^~ ww - ^2^(2.-1) R -1 



MATCHING * o * 



JiV- 



* and ^ , ^f tc ' l (A, x) are matching by construction. The monodromy of * at A = 0,x coincides with 
that of ^fj^ tch {\,x). In order to compute it, we write: 



*o(M) = A 1 2 (A* ~ !) ? $o(a*), $o(a*) 



<Pl ¥>2 

6 6 



$ satisfies (49) with u = 1. It is expressed in terms of two independent solutions <p 2 of the 
Gauss hyper-geometric equation (see Appendix 1): 

Ml - A*) ^ + (1 + c - (a + [b + 1] + 1) p) ^ - a(b + 1) <p = 0, 



where, 



From ip i and y>2 we compute: 

6-1 

r 



1,2. 



We need a compete set of solutions at p = 0, 1, oo. Since a— 6 is an integer, we are in a logarithmic 
case. We briefly explain some preliminary facts. Let us consider a Gauss Hyper-geometric equation 
in standard form: 

p(l-p)^+[j-(a + /3 + l) p]^--af3<p = 

(a, /3,7 here are not the coefficients of (PVI)! We are just using the same symbols only here). 
Logarithmic solutions at p — may occur only if 7 6 Z. Logarithmic solutions at p = 1 may occur 
only if a + /3 — 7 e Z. Logarithmic solutions at /U = 00 may occur only if a — (3 G Z. Several sub-case 
must be distinguished, and this is not the place to discuss them. 

In our case a = a, ,8 = 6 + 1 and 7 = c + 1. Therefore a — j3 — —2. This case is logarithmic. 
Two independent solutions are: 

„(°°) _ ..-£„ ( a a „, 1 i 1 1 a „,. M _ r ;„ „™ 1 ,.-«- 2 , ^ — - - 1 



(pi 00 ' = p p gi ( p,p - 7 + 1, 1 + fi - a; - J = [ in our case ] p gi I /3, /3 - 7 + 1, 3; - J , 
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^ =fi- p F 7 + 1,1 + -a; ij = [ in our case ] pT^F ^,^-7 + 1,3; ij . 

Here F is the Gauss hyper-geometric function and 171 is a logarithmic solution introduced in Nor- 
lund's paper [23], page 7: 

gi (u,v,w;z) = yV-l)"- 1 ^ - 1)! 1 + F(u,v,w; z) ln(-z) + 

+ £ [V(l - u - n) + V(« + n) - V(™ + n) - V(l + «)] « n , 

„=o n 'W« 

where 

^(z) := -f-lnT(a;), a; ^ 0, -1, -2, -3, ... ; Id < 1, ln(-z) < for - 1 < z < 0. 
ax 

Note that the first sum J2n=i(~ 1)™ _1 (^ — 1)! ^"^ u ^~" is a polynomial in 1/z. 
If #0, 0^ ^ Z, we fall in the non-logarithmic cases at fi — 0, 1. Thus: 

(o) ^ f 0q X 1 #0 . 6x , 3 \ 

^ - F ( V T + T-2'T + T + 2' 1 + ^' 

(o) -Bar? ( Oq x 1 0o . 0s . 3 
(i) t-, / 0o , 0s 1 0o . 0k . 3 \ 

w =F U + ^"2'T + T + 2' 1 + ^ ;1 - M J' 

(1) Ci \-8^Tp{ 0o 0x 1 0o 02 . 3 \ 

^ ; = (i-m) ^( v y-y-2'y-T + 2' 1 -^ ;1 " /i J- 

(00) _£o_»x_3 { d x 3 O X 3 1\ 
^ 2 2 2 51 U + ^ + 2'~T + T + 2' 3; mJ' 

( TC)= F C«o + ^ + 3 00 0s 3 n 

y2 ' V 2 2 2 2 2 2 V/ 

The connection matrix between [y^, ip^] and [^i^, ¥>2^] i s a standard one: 
[^ 0) ,^ 0) ] = [^\^]C 0U < arg (1 -n) < 27r; 

where Cqi is (36). 

On the other hand, the connection matrix between [ip^\ ip^} and [ip^, ip^} is computed in 
[23]. Our specific case falls in the case a — (3 = —p, p > integer. From [23], page 27, case 11, 
formulae (1) and (2), we deduce the connection formulae: 

i-„r-"F (0,0- 7 + 1. 1 + 0- a; I) = + 



(-rt° _s Ji (/?,0-7 + 1,1 + /?-«; i) 
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Here |arg (— fj,)\ < ir, ln(— /j.) > for — oo < /j, < — 1. The branch cut in the /z-plane is [0, +oo). 
the minus signs, we choose — fi = e~ 47r ^t, then < arg^t < 2n. Moreover, in our case (— 1) Q ~ 
With this preparation, we can write the connection matrix: 



[^ oo) ,^ oo) ] = [^ 0) ,^ 0) ] Cooo, < ar gAt < 27r; 



where Cqoo is (35) 
In order to wi 



In order to write ^f^ tch 

(oo) 

Aoo) 
S2 



(oo) 
L 

oo) 
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[(0 O + e x y - i][(e - e x y - 1] 



i + o ( - 



+ 2 In /x 



M 2 — 1) 2 ¥?2 



oo) 



1 + 



As for 
= 1. 



in terms of the tp\°°\ let us compute the behavior for fi — > oo of 
. We have: 



1 + - 



M 2 - 1) 2 4i = A* 2 hi^ 



M 2 (A» " 1) 2 * Q 



A* 



1 + - 

A* 



l + O 



+ o( M - 2 ). 



Therefore: 



/ (oo) (oo) 



/ o( - 

A* 



[J? 



i o 

i? In/x 1 



We conclude that the matrix *o used in the ^oc/t ^ */jv matching is: 

C, C := 



2 



(oo) (oo) 

*o = A^>-l)^ ( ^oo) %) 

, SI 



-J? § 
5 
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As we have already remarked, we do the transformation ij/^* c;i i_> ~$>^ tch CiN , to compute 
the x-independent monodromy at A = 0, x. This means that we have to compute the monodromy 
of: 

Kn^Cin- 1 = K (x)^ Q (A/x) = 

«0 



(oo) (oo) 

1 V 

^ e 2 



xj \x J V4T £ 



>2 



x / \ x 



2 



In order to do this, we observe that from the definition of ipf\ip^ if follows that: 

for A — > 0. From the definition of ip^\ip^\ it follows that: 

Ko{x) (i)' Q" 1 ) 2 (tf> ^^^(^(^^(A-^KA-x)^, 
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for A — > x. It is not necessary here to write explicitly the invertible matrices tpQ N (x) and ipi Nx - The 
above construction implies that ^fj^^CiN^ 1 has monodromy matrices at A = 0,x respectively 
given by: 

M = C~ 1 (C oo0 ) exp{i7r6> o-3 

M x = C-^C^o)" 1 ^!)" 1 exp{iw6 x a 3 } CiC^C 

These matrices coincide with the monodromy matrices of \P Cin 1 - 

As a last simplification, we consider the transformation M i->- CMC -1 . We obtain the result of 
theorem 3, case c). Namely: 

M, = (Coco) 1 exp{i7r6» cr 3 } Coo , M M = ^ _ ^ , 

M i = (^27ria l)' M x = (C oo0 ) 1 (C 01 ) 1 exp{i 7 r6' 2; cr3} CoiCooo- 
We observe that r does not appear in the monodromy matrices. 

8 Appendix 1 

Proposition 2 Let Bo, B\ fee 2 x 2 matrices such that 

Eigenvalues (Bo) = 0, — c, Eigenvalues (B\) = 0, c — a — b. 
and Bo + B\ is either diagonalizable: 

Bo + B\ = ^ ^ (it may happen that a = b), 

or it is a Jordan form: 

Bo + B 1= ( - a 1 



-a)' 

Then, B and B\ can be computed as in the following cases. Let r, s be any complex numbers. 

1) Diagonalizable case. 

Case a ^ b: 

/ a(b-c) r 



a) 
b 


— r 


))21 


b(b-c) 
a — b 




-r N 


(0 


c-b, 


(c- 


a — r 


\ 






Bo := ab(a a S c )( c -b) bjc-a) . B 1 = ( «"» h(h _ r , ] , r + (98) 

\ r(a-b) 2 a-b ) 

ifa = 0: S =(o ^ ) , B 1 = ( ^ ^'J. (99) 

//6 = 0: £o=(- C J), Bi=^ Q - 'j. (100) 

Jf« = c*6: 50^(7 j) , = (j ^ . (101) 

//6 = c^a: S =(° Bl =(- 7). (102) 

Cases (99)-(102) are already included in (98). 
Case a — b. We have two sub-cases: 

Ifa = b = c: B = ^ f (e+ j) g i = ( s ( c + s ) _ c ^ s ) ' 1!}:n 

rf„ = b = (): Up = ( "f (c +f) Bi = -B . (104) 
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The transpose matrices of all the above cases are also possible. 

2) Jordan case. 

For a ^ and a ^ c we have: 



r(r+c) \ / _ 1 _ r(r+c) 

Bo = ( , N ,51=, 

, a(c — aj — c — r / \a(a — c) c — a + r 



For a = 0, or a = c, we have two possibilities: 

B = ( n ^ I > 



or 



Bn = 



— c r 




Si = 



—a 1 — r 
-a + c 

c — a 1 — r 



Proposition 3 iei Bp and B\ be as in Proposition 2. The linear system: 



d_ ( if 



Bo , Bi 



z - 1 



may be reduced to a Gauss hyper-geometric equation, in the following cases. 
Diagonalizable case (i.e. from (98) to (104)): 

z(l-z)^ + (l + c-(a+[b+l] + l) z) ^-o(6+l) if = 0. 



dz 



(105) 

(106) 
(107) 



(108) 



The component!; is obtained by the following equalities, according to the different cases of Proposition 
2. 



Cases (98) (99) (100) (101) (102): 



dip ( b-c 
Z{l - Z) TZ a Z + * 



(109) 



Case (103): 



Case (104): 



/-, x dtp , . 
z{l — z) — + (c + s - C Z) ip 
dz 



z{l-z) ^ + (c + s) tp 
dz 



Jordan case (105): The equation for ip is in Gauss hypergeometic form only when r — —a. In this 
case, the matrices (105) are: 



The equation is: 



^ a(c-a) a — c) ' ^ \a(a — c) c — 2a /' 



z(l - z) ^ + [1 + c - (2a + 1) z] J - a 2 <p = 0, 



t, = z — + a ip. 

dz 



Jordan case (106): for r — 1 we get: 
For a = 0, 



«( 1 -^ + ( 1 + c--)^ = () - 



dz 
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Lp = 0. 



Forc = a, z(l - z)§ + (1 + a - (2a + l)z) J + (-a 2 - ^ 

T/ie Gauss form appears only when a = 0. 
Jordan case (107): for r = 1 we get: 

. d 2 tp , . dw c 

Fora = 0, ,(1 - *)-^ + (1 + c - *) J + — „ = 0. 

. d 2 LD , , . . aV n 

For c = a, z(l - z)-r-£ + (1 + a - (2a + l)z)-r- - a 2 <^ = 0. 

dz 1 dz 

For r ^ 1, we don't get a Gauss hyper- geometric form for the equation of ip in both cases (106) and 
(107). Nevertheless, the matrices are in upper triangular form, so the equation for £ is solvable by 
elementary integration. 

Jordan case (105): The equation for £ is in Gauss hypergeometic form: 

Z{Z ~ + + C ~ 2(a + 1)Z) ^ ~ a{a + 1)e = °' 



dz 



ip(z) 



a(a — c) 



z(z — 1)- — h (az — c — r)£ 



As a compendium to the above proposition, we recall that any irreducible representation of 
7r 1 (CP 1 \{0, 1, oo}) M> GL(2, C) can be realized as the monodromy group of a Riemann (or Gauss) 
equation. A reducible representation (namely, Mq,Mi,M 00 are in upper triangular form) can be 
realized by the monodromy of a 2 x 2 Fuchsian system 



dY 

dz 



Bp Bi 

z z — 1 



where Bp, B\ are 2x2 upper triangular matrices. We also state the following: 
Lemma 12 Consider a 2 x 2 linear system: 



dY(z) 
dz 



A(z) Y, A(z) 



a(z) b(z) 
c(z) 



such that A(z) is meromorphic, with poles a\, a 2 , on, oo. The monodromy group is generated 
by N upper triangular monodromy matrices: 



Mi 



A« J**)' 
A« )' , " 1 ' 2 '-' 



w/iere A^A^-R^ are constants (i.e. they are independent of z) given by: 

A^ = exp{27ri Res a(z)| ai }, A^ = exp{27ri Res c(z)| ai }, 



Z CXp{27Tz} 



ds b(s) 



U2(S) 

m(s) 



ds o(s) 2 [ | . 
«i(«) 



i?W depends on a chosen non-singular base point Zq, but not on z. One of the non zero i?W caw fr e 
pwf egita^ to 7 7 fry re-defining Y M> FC, C = ^ ^ " 

Proof: Let us write F = ^ ^ ^ and the equation in the form: 



dy 



".'/i , • , dy 2 
— = a(z) yi + b{z)y 2 , = c(z)y 2 . 
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Let zq 7^ di,oo, i = l,...,iV. The second equation has solution y2{z) = C 2 w 2 , where U2 := 
exp{ c(s)}, C2 G C. The first equation becomes: 

^=a(z) yi + C 2 b{z)u 2 {z). (110) 

We solve the first equation by variation of parameters. Let u\(z) — exp{ f* ds a(s)} be a fundamental 
solution of the homogeneous equation = a{z)y\. We look for a solution of (110) of the form 
2/2(2) = w(z) ui(z). Substitution gives: 

— = C 2 b{z) — ^ w(z) = C 2 v(z) + Ci, Ci G C. 
az Ui(z) 

where: 

«(*)= r^^rS 

The general solution of (110) is yi(z) — C\ u\(z) + C2 v(z) u\(z). Then, a fundamental solution 
for the initial system can be chosen to be: 



Y(z) 



'ui(z) v(z) ui(zY 
u 2 (z) 

We compute the monodromy for (z — aj) h-> [z — ai)e 2m . We have: 

u\(z) i->- A^ ui(z), xf = cxp{27U Res a(z)\ ai }, 

u 2 {z) h-> xf Ul (z), xf = cxp{27ri Res c(z)\ ai }. (Ill) 
By linearity, the vector solution f v ^ "i( z ) j j s transformed into a linear combination of two 



u 2 (z) 

independent vector solutions: i?W ^^q^^ + (z|^)' ^ oreover ' must coincide 

with xf, because - by virtue of (111) - u 2 (z) n> A^u 2 (z). Namely: 



v(z) m(z) 
u 2 (z) 



Thus: 



'ui{z) v(z)u 1 (z)\ (xfu x {z) Xf v(z) Ul (z) + i?« Ul {z)' 

r, '- ,= 1 ' 

u 2 (z) J V 4° "2(2) 



Y(z) 



Xf JJ« 
A« 



Let Ci be a small loop around a.j. To find let us observe that u\(z) v(z) i-> A^ ui(z) v(z), 
where: {;(*) = + ATi(z), := /; oxp ^^> ^ 6(a) ^g}. Thus: 

/ Ul (z) /AWui(z) A< 4) «(z) Ul (z) + Jfi(z)ui(«)' 

y(z) = ^ 

V u 2 (z) / V A 2 * } u 2 (z) 

We must have: xfv(z)ui(z) + Ki(z)u±(z) = xfv(z)u±(z) + R^ui(z); namely: 

+ (Xf - Xf)v(z). 

□ 
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9 Appendix 2: Formal Asymptotic Expansion 



1) We consider systems (62): 



dY 

dz 



ft + — + V — 

n=2 



Y :=D(z)Y, ft = diag(wi,w 2 , w n ), 



with 7^ Wj, for i 7^ j. We introduce a gauge transformation Y = G(z)Y, such that: 



Y 

dz 



G-\z)D{z)G{z)-G-\z) d ^- 



dz 



Y, 



be in the simple form: 



dY 

dz 



ft 



fti 



Y, ft, fti diagonal. 



Formally, we write G{z) as: 



G(z)= 7 + — + % + ..., z^oc. 

z z 

If the G n 's can be determined, we get the formal solution: 



Y(, 



oc „ 



exp {z ft + fti m(z)} , z — > 00 



For a sector of angular with ir + e, e > sufficiently small (but finite and non zero), there exists a 
unique solution Y(z) with the above asymptotic expansion [3]. 

In order to determine G n and fti, we solve D(z)G(z) — d z G(z) = G(z)(fto + ftiz -1 ): 



K n=l / \ 71=1 / n=2 



(n - l)G„_l 



We identify equal powers of z . From the power 1/z we get: 



n 1= diag^,...,^), := (ni)« = (£>i)«, (Gi)y = — 



From the power l/z 2 , for « 7^ j, we compute (G2)ij, and for i = j we compute: 



(Gi)ii = -(D 2 )u - y^(£ ) i)ifc(Gi)fcz- 



From the power 1/z" we get: 



(G„_i)jj — - ^ < — (£>„ + D n _iGi + ... + D 2 G n -2) ii — ^(£ , i)ife(G„_i)fei > , 



(G n )ij 



u>i — 00 i 



(!)_,,« 



- {D n + D„_iGi + ... + D 2 Gn-2)«}, * ^ j 



2) We consider the system (63): 



dY 

dz 



z ^ z n 

n=2 



Y :=E(z) Y, A= diag (Ai, A 2 , A„). 
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with Xi =/= Xj for i =/= j. We introduce a gauge transformation: 



in order to reach the simple form: 

dY 

dz 



Y(z)=K(z) Y(z) 



K- 1 {z)E(z)K(z) - K-^z)^^- 

dz 



Y 



2 a a A l 

x z A z + xA H 



Y, Ai diagonal. 



Formally, we write: 



n=l 



Provided that we can determine the matrices K n , we obtain the formal solution: 



Y(z) 



oo 

i+y ^ 



n=l 



cxp ^ —A z + xA z Ai In a; 



For a sector of angular with ^ + e, e > sufficiently small (but finite and non zero), there exists a 
unique solution Y(z) with the above asymptotic expansion [3]. 

In order to determine K n and Ai, we solve E(z)K(z) - d z K(z) = K(z)(x 2 Az 2 + Az + Sliz -1 ): 



(n - l)K n _i 



71=1 



= (' + g^)(^ + A= + ^)- 



We identify equal powers of z^ 1 . 
Power z. It is an identity: x 2 A = x 2 A. 
Power z . 

xA + x 2 AK\ = xA + x 2 K 1 A 
This means that K\ is a diagonal matrix. 
Power 1/z. We obtain the equation: 



[A, Ki] = 



x 2 [A,K 2 }=A 1 -E 1 



Power 1/z 2 . We obtain the equation: 



(Ai)ii — (Ei) a, (G 2 )ij 



( E l)ij 



x 2 (Xi - Xj) 



x 2 [A, K 3 ] = x[K 2 ,A] + KiAi - E X K X - K x - E 2 = 



Ei-A l 



+ KiAi - EiKi - Ki - E 2 



Thus: 



(K 3 )i 



1 



x 2 (Xi — Xj ) 
Power 1/z 2 . We obtain the equation: 



(Ki) u = -(E 2 )i 



(Ei)ij - + (E 2 )jj - (E 2 ) 



x 2 [A, Ki] = x[K 3 , A] + K 2 Ai - E X K 2 - E 2 K 1 - 2K 2 - E 3 



The diagonal part gives: 



2(K 2 )ij = [(E 2 ) u } 2 - (E 3 )u-Y,(Ei)ik(K 2 ) k i. 

k=£i 



The non diagonal part gives (K^)^, i ^ j. 

We content ourselves with these results, namely the determination of K\ and K 2 . With the same 
procedure, we can determine all the K^s. 
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10 Appendix 3: Birational Transformations 

All the solutions of (PVI) of the form: 

y{x) = b + hx + b 2 x 2 + ... + b N x N + ... (112) 

are obtained from the matching procedure of sections 5,6,7. By this we mean that solutions of type 
(112) are the solutions given by the matching procedure, or they can be obtained from solutions given 
by the matching procedure via one of the birational transformations of [24] and the transformation 
(22). 

Birational transformations are symmetries of (PVI), namely invertible transformations: 

^ ) = iSpi' x '=wy <*■«■■•■•«-) 

such that y(x) satisfies (PVI) with coefficients 0o,8 x ,0\,0oo an d variable x, if and only if y'(x') 
satisfies (PVI) with coefficients 8' Q ,8 X ,8' 1 ,8' 00 and variable x' . The functions P, Q are polynomials; 
p,q are linear; the transformation of the # M 's is an element of a linear representation of one of 
the following groups. Permutation group; the Weyl group of the root system D 4 ; the group of 
translations v := (v\, v 2 , v 3 , V4) M> v + ej, j=l,2,3,4 (where e\ — (1, 0, 0, 0), e 4 — (0, 0, 0, 1)). 7 

* Permutation group: 

x 1 : 0i = o , o' Q = e i; e' x = e x , e' 00 = e 00 ; y '{x) = i- y (x), x = i-x'. 



„2 



1 1 



X 



1, 6^ = 60 + 1; 0[ = 0i, e' x = x; y '( x ) = —, x 



x 3 : e' = e x , e' x = e ; 0[ = 1 , 0^ = 0^, y'(x) = --^-, x= r 

x — 1 



* Weyl Group: 

w 2 : 



8[ — -81; 9' — 6 Q , 6' x 



% + 81 + X + #oo 8 + 81 — 8 X — 8oo 
2 h ° 1 = 2 + 1 ' 



a , "0 - <?1 + Vx - 80c . , a i &0 - VI - Ox + t/oo . 1 
U x - 2 h ' 00 = 2 



w 3 : 



W4: 



O'oo — 2 — ^oo5 0'a — 00, 0'x — u x, 0'\ — 0\. 
J'oo = 2 — 0oo; ^ = 2 — $2,; = # 0j #'1 = 



The variable x' — x, but y'(x) is quite complicated and will not be given here (see [24]). 
* Shift lj : v t-> v + ej : 



h ■■ 


O'o 


= 0o + l, 


8[ = 8 1 + 1; 


0'x - 0x, 


0'oc — 0oc 


h ■■ 


O'o 


= 0o + h 


8[ =0i-l; 


0' x = 0x, 


0'oC = 0oG 


h: 


0'x 


= 0x + l, 


o'oc = + 1; 


0'q = 00 


, 8[ = 81 


u ■■ 


0'x 


= 0x + h 


0'co = Occ - 1; 


6' = 8 . 


, 8[ = 81 



The variable x' — x, but y'(x) is quite complicated and will not be given here (see [24]). 



7 We note that 0i , X , 0q are defined up to the sign. Moreover, 6qc is defined up to #00 i— >■ 2 — $oo, and for this 
reason symmetries arc described in [24] in terms of Xoc := #oc — 1 . 
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For the Taylor solutions, we have a = 9\ — 8^. Denote a' := 8[ — 9'^. If we start from 
a Taylor solution constructed in sections 5, 6, 7 by means of the matching procedure, developed 
for a = 9\ — #00, then the birational transformations allow to obtain the solutions defined for 
a = ±(#i ± \ T $oo) + n, n G Z. This is a consequence of the following actions: 

li and I4 : a i-> a' = a + 1. I2 and I3 : a ^ a' = a — 1 

W3 ■ h ■ h : o 1 y (7 1 == 9\ -\- #00. u>i ■ W3 ■ h ■ h : a a' = —a. 
Note also that (PVI) is invariant for 8\ n> — B\. This maps a i->- — cr. Other actions are: 

wi : a' = -(81 +#oo); w 2 : <?' = cr; w 3 : cr' = #1 +6>oo -2; w 4 : cr' = 6>i + 0,^ - 2. 

x 1 : cr' = # -#oo; a; 2 : cr' = 6>i - O ; a; 3 : cr' = #i - 6^ = a. 



11 Appendix 4: Examples of Taylor Expansions 

We give the solutions of (PVI) of the form y(x) = X^^Ln ^" ^° ^ depending on the value of 
the coefficients 9q,9\,9 x , 9^. Similar examples can be constructed for solutions y(x) = Y^Li bn x n , 
bi ^ 0, by the symmetry (22). 

We observe that, in general, the free parameter appears starting from some power x m . If we 
truncate the series at a; m_1 , we cannot see it. 

We always denote b n the coefficients, though they change case by case. 

* Example 1: (PVI) always has one solution (18) when 8^ — Q\ Z, and one solution (21) 
when 9 00 + 9 1 ^ Z. 

* Example 2: Case 9 1 + 9 X = 0. 

i) There are solutions (19), if 9 X — ±9 0} Boo 1- 

ii) There is a solution (18), defined for 9^ ^ 1, ^±1, where n G Z. The condition B x ^ 2«±! 
follows from the condition 6*i — #oc ^ n in (18), when 8\ = —8^. The denominators of the coefficients 
b n vanish for half- integer 8^, and for 9^ = 1. 

iii) For 9^ — 2 " 2 t " 1 , solutions y(x) = X^Lo ^n 21 ™ m ii) are n °t defined. On the other hand, solutions 
(19) are defined, provided that 8 X = ±#o- Moreover, for a given = , there may be solutions 
equivalent to (19), provided that 8 x ±8q assumes some integer values. For example, consider #oc = §• 
We have (19): 

/ 3 1 \ 00 3 

y(x) = -2 +ax + (0 2 - 1 + - a - - a 3 J x 2 + ^ b n {a; B , -) x n , 8 X = ±8 Q ; 

~ ' n=3 

and a solution equivalent to (19): 

00 

y{x) = 4 - (2 ± 8 ) x +ax 2 + ^ b n (a; 9 ) x n , 9 X = ±(1 + B ), ± (1 - 8 a ). 

n—3 



iv) If = 0, 1, we do not have any solution of the desired form, except for the singular solutions 
y = 0, 1. If 800 — 2 we have a solution equivalent to (20). The parameter a is the coefficient of a; 3 
(coefficients of x°, x and x 2 have no parameter). 

* Example 3: Case 8^-81 = 0. 

i) We have solutions (19) 

ii) We have the solution (21), with the substitution: 9\ = 9^. This is defined for 9^ ^ 1, 2 " 2 +1 . 

iii) For = 2 " 2 +1 we find solutions of exactly the form (19), when 9 X = ±80. Moreover, for any 
given 800 = 2 " 2 t " 1 , there may be solutions equivalent to (19), provided that 9 X ±80 has some integer 
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value. For example, consider — —\- There are solutions (19): 

2 ^ 1 

V{x) = - +ax + 2_^b n (a;6 ,--) x n , X = ±9 ; 



and the equivalent solutions: 

V{x) = \~ x +ax2 + £ b ^ a ' oK, Ox = ±(1 + Oo) or ± (1 - O ). 



n—3 



Another example: consider 9^ = |. In this case we just have (19), or the singular solution y=l. 

iv) For #oo = 0, 1 we don't have solutions of the desired form, except for the singular solutions 
y = 0, 1. If 6*oo = 2 we have a solution equivalent to (20). The parameter a is the coefficient of x 3 
(coefficients of x , x and x 2 have no parameter). 

* Example 4: Case 61 — 600 = — 1- 

i) We have the solution (21), defined for 0^ 7^ 1, where n G Z. 

ii) For any #oo = 2 1 " 1 , solutions (19) are defined, when 9o ±9 X = 0. Moreover, solutions may exist 
equivalent to (19) by symmetry, provided that #0 ± 9 X is some integer. 

iii) For 6^ = 1 we are exactly in the case (20). 

iv) For 9oo — 0, 2 we have solutions equivalent to (20): 

y{x) = 2 x +a x2 + ^-e x f- m e^e x f-m-ol) + a fl§-*, + 3j x3 + 

* Example 5: Case #i — ^ = 1. 

i) We have solution (21), with the substitution of 6\ = 9^ + 1, and defined for 9^ ^ 1, 2 " 2 +1 , n € Z. 

ii) Equivalent to (19) by symmetry, we have the solution: 

/ \ 2 (#00 + 1)(2 ± O ) . 2 , 1 1 a Q \ n 
y( x ) = - 1 2i — + ^Ta T\ x +ax z + yb n (a,e oo ,0 o )x n , 

1 — Poo O^foo — 1J 

floo ^0,1; X = ±(1-0 O ) or ±(1 + O ) 

The two signs in the coefficient of x depend on the choice 9 X = ±(1 — 9 ) or ±(1 + O ) respectively 
Similar change of signs occur in all the coefficients b n . 

iii) If 6>oo = 0, the solution in case ii) is not defined (denominators in the coefficients b n vanish). We 
have anyway a solution equivalent to (20): 

02 _ 02 _ o 00 

y(x) = 2 x +ax 2 + V b n (a; 0o,0oo)z n . 

6 ^3 

If 6*oo = 2, we have a solution equivalent to (20). The parameter a is the coefficient of x 4 (no 
parameter in lower powers of x). 

iv) If 6>oo = 2 "+ 1 the solution in i) is not defined. Solutions exist equivalent to (19), provided that 
9 X ± 6*o is some integer. 

v) For 6*oo = — 1 we have a solution equivalent to (20): 



y(x) = l + ax 2 + a{6x °° +3 K 3 +Y / b n (a;9 ,8 x )x". 

n=4 

vi) For 0oo = 1, solutions of the desired form do not exist, except for y = 1. 
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We could proceed at our pleasure, choosing any value 9^ ± Q\ integer. We would always have 
solutions of three kinds. 1) One out of the two solutions (18) and (21). 2) Solutions equivalent to 
(19) - at least when is half integer - provide that X ± do is some integer. 3) Solutions equivalent 
to (20), for 9oo equal to some integer. 
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